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Abstract 

A NNLO analysis of certain logarithmic expansions, developed for precision stud¬ 
ies of the evolution of the QCD parton distributions (pdf) at the Large Hadron 
Collider, is presented. We elaborate on their relations to all the solutions of the 
DGLAP eqnations that have been hitherto obtained from Mellin space, to which 
are equivalent. Exact expansions, equivalent to exact solutions of the equations, are 
constructed in the non-singlet sector. The algorithmic features of our approach are 
also emphasized, since this method allows to obtain numerical solutions of the evo¬ 
lution equations with the same accuracy of other methods, based on Mellin space, 
and of brute force methods, which solve the equations by finite differences. The 
implementation of our analysis allows to compare with existing benchmarks for the 
evolution of the pdf’s, useful for applications at the LHC, and to extend them sig¬ 
nificantly in a systematic fashion, especially when solutions that retain logarithmic 
corrections only of a certain accnracy are searched for. 
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1 Introduction 


Precision studies of some hadronic processes in the perturbative regime are going to be 
very important in order to confirm the validity of the mechanism of mass generation 
in the Standard Model at the new collider, the LHC. This program involves a rather 
complex analysis of the QCD background, with the corresponding radiative corrections 
taken into account to higher orders. Studies of these corrections for specific processes have 
been performed by various groups, at a level of accuracy which has reached the next-to- 
next-to-leading order (NNLO) in the QCD coupling constant. The quantification of 
the impact of these corrections requires the determination of the hard scattering of the 
partonic cross sections up to order with the matrix of the anomalous dimensions of 
the DGLAP kernels determined at the same perturbative order. 

The study of the evolution of the parton distributions is critical for the success of this 
program and may include both a NNLO analysis and, possibly, a resummation of the 
large logs which may appear in certain kinematic regions of specific processes p. The 
questions that we address in our work concern the types of approximations which are 
involved when we try to solve the DGLAP equations to higher orders and the differences 
among the various methods proposed for their solution. 

The clarification of these issues is important, since a chosen method has a direct im¬ 
pact on the structure of the evolution codes and on their phenomenological predictions. 
We address these questions by going over a discussions of these methods and, in partic¬ 
ular, we compare those based in Mellin space and the analogous ones based in x-space. 
Mellin methods have been the most popular and have been implemented up to NLO and, 
very recently, also at NNLO j2j. We remark that x-space methods based on logarithmic 
expansions have never been thoroughly justified in the previous literature even at NLO, 
in the case of the QCD parton distributions (pdf’s) jSlIllEI- 

We fill this gap and present exact proofs of the equivalence of these methods - in the 
case of the evolution of the QCD pdf’s - extending a proof which had been outlined by 
Rossi |fi] at LO and by Da Luz Vieira and Storrow |3| at NLO in their study of the parton 
distributions of the photon. 

In more recent times, these studies on the pdf’s of the photon have triggered similar 
studies also for the QCD pdf’s. The result of these efforts was the proposal of new 
expansions for the quark and gluon parton distributions uw which had to capture the 
logarithmic behaviour of the solution up to NLO. Evidence of the consistency of the ansatz 
was, in part, based on a comparative study of the generic structures of the logarithms that 
appear in the solution using Mellin moments, since the same logarithms of the coupling 
could be reobtained by recursion relations. 

In this work we are going to clarify - using the exact solutions of the corresponding 
recursion relations - the role of these previous expansions and present their generaliza- 
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tions. In particular we will show that they can be extended to retain higher logarithmic 
corrections and how they can be made exact. It is shown that by a suitable extension of 
this analysis, all that has been known so far in moment space can be reobtained directly 
from x-space. In the photon case the Da Luz Vieira-Storrow solution |2] can be now 
understood simply as a first truncated ansatz of the general truncated solutions that we 
analize. However, the analysis presented here is limited to the QCD pdf’s, while a similar 
study of the pdf’s of the photon will be discussed by us in a forthcoming paper. 

We introduce the notion of first, second, and so on, truncated solutions of the DGLAP, 
that correct the LO behaviour, from which the nature of the expansions and the contri¬ 
butions retained by the various approximations will appear quite clearly. 

Our work is organized as follows. After defining our conventions, we bring in a simple 
example that shows how a non-singlet LO solution of the DGLAP is obtained by an x- 
space ansatz. Then we move to NLO and introduce the notion of truncated logarithmic 
solutions at this order, moving afterwards to define exact recursive solutions from x-space. 
In all these cases, we show that these solutions contain exactly the same information of 
those obtained in Mellin space, to which they turn out to be equivalent. The same 
analysis is then extended to NNLO. Our approach lays at the foundation of a numerical 
method -based on x-space - that solves the NNLO DGLAP with great accuracy down 
to very small-x (10“^). The method, therefore, not only does not suffer from the usual 
well known inaccuracy of x-space based approaches at small-x values jHj, but is, from our 
viewpoint, a complementary way to look at the evolution of the pdf’s in an extremely 
simple fashion. We conclude with some comments concerning the timely issue of defining 
benchmarks for the evolution of the pdf’s, obtained by comparing solutions extracted by 
Mellin methods against those derived from our approach, in particular for those solutions 
which retain accuracy of a given order in (0 (q!s) accurate solutions), relevant for precise 
determination of certain NNLO observables at the LHC. 


2 Definitions and conventions 


Before we start our analysis it is convenient to define here our notations and conventions 
that we will use in the rest of the paper. 

We introduce the 3-loop evolution of the coupling via its /3-function 
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and its three-loop expansion is 
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are the coefficients of the beta function. In particular, P 2 pum and /ds m are in the 
MS scheme. We have set 

7Y2 _ 4 

Nc = 3, Cp = = -, Tf = TpUf = -Hf, (4) 

where Nc is the number of colors, rif is the number of active flavors, that is flxed by 
the number of quarks with niq < Q. One can obtain either an exact or an accurate 
(truncated) solution of this equation. An exact solution includes higher order effects in 
tts, while a truncated solution retains contributions only up to a given (flxed) order in a 
certain expansion parameter. The structure of the NLO exact solution of the RGE for 
the coupling is well known and relates ctsil^i) in terms of as(/i 2 ) via an implicit solution 
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where = Q!s(/r^)/(47r). The truncated solution is obtained by expanding up to a 

given order in a small variable 
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where the dependence is shifted into the factor L = ln(/i^//r|), and we have used 
a /3-function expanded up to NLO, involving /3o and /3i. Exact solutions of the RGE 
for the running coupling are not available (analytically) beyond NLO, while they can 
be obtained numerically. Truncated solutions instead can be obtained quite easily, for 
instance expanding in terms of the logarithm of a specific scale (A) 
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and where is calculated using the known value of as{mz) and imposing the 

continuity of at the thresholds identified by the quark masses. 
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3 General Issues 


For an integro-differential equation of DGLAP type, which is defined in a perturbative 
fashion, the kernel P{x) is known perturbatively up to the first few orders in approx¬ 
imations which are commonly known as LO, NLO, NNLO. 

The equation is of the form 

= P{x,Q^)^ f{x,Q^), (9) 

with ^ 

a{x)®b{x)= [ —a{y)h{x/y), (10) 

Jo y 

and the expansion of the kernel at LO, for instance, is given by 

P(x, Q^f° = (^^) ( 11 ) 

In the case of QCD one equation is scalar, termed non-singlet, the other equation involves 
2-by-2 matrices, the singlet. In other cases, for instance in supersymmetric QCD, both 
the singlet and the non-singlet equations have a matrix structure na. Except for the LO 
case, exact analytic solutions of the singlet equations are not known. However, various 
methods are available in order to obtain a numerical solution with a good accuracy. These 
methods are of two types: brute force approaches based in x-space and those based on 
the inversion of the Mellin moments. 

A brute force method involves a numerical solution of the PDE based on finite differ¬ 
ences schemes. One can easily find a stability scheme in which the differential operator 
on the left-hand-side of the equation gets replaced by its finite difference expression. This 
method has the advantage that it allows to obtain the so called “exact” solution of the 
equation at a given order (LO, NLO, NNLO). The only approximation involved in this 
numerical solution comes from the perturbative expansion of the kernels. Solutions of 
this type are not accurate to the working order of the expansion of the kernels, since they 
retain higher order terms in ^ A short-come of brute force methods is the lack of an 
ansatz for the solution, which could instead be quite useful in order to understand the role 
of the retained perturbative logarithms. The use of Mellin inversion allows to extract, in 
the non singlet case, the exact solution quite immediately up to NNLO. 

The Mellin moments are defined as 

^It is common, however, to refer to these solutions as to the “exact” ones, though they have no better 
status than the accurate (truncated) ones. In principle, large cancelations between contributions of higher 
order in the perturbative expansion of the kernels beyond NNLO, which are not available, and the known 
contributions, could take place at higher orders and this possibility remains unaccounted for in these 
“exact” solution. The term “exact”, though being a misnomer, is however wide spread in the context of 
perturbative applications and for this reason we will use it throughout our work. 
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a{N) = / a{x)x^ 

Jo 

and the basic advantage of working in moment space is to rednce the convolution product 
0 into an ordinary product. For instance, at leading order we obtain the LO DGLAP 
equation © in moment space 


df{N,a,) 


da,. 
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which is solved by 
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where we have used the notation a = a(Q‘^) and ckq = «(Qo)- ^^is point, to contruct the 
solution in x-space, we need to perform a numerical inversion of the moments, following a 
contour in the complex plane. This method is widely used in the numerical construction 
of the solutions and various optimization of this technique have been proposed [21]. We 
will show below how one can find solutions of any desired accuracy by using a set of 
recursion relations without the need of using numerical inversion of the Mellin moments. 


3.1 The logarithmic ansatz in LO 


To illustrate how the logarithmic expansion works and why it can reproduce the same 
solutions obtained from moment space, it is convenient, for simplicity, to work at LO. 
We try, in the ansatz, to organize the logarithmic behaviour of the solution in terms of 
as and its logarithmic powers, times some scale-invariant functions An(x), which depend 
only on Bjorken x. As we are going to see, this re-arrangement of the scale dependent 
terms is rather general for evolution equations in QCD. The number of the scale-invariant 
functions An is actually infinite, and they are obtained recursively from a given initial 
condition. 


The expansion that summarizes the logarithmic behaviour of the solution at LO is 
chosen of the form 
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To determine An{x) for every n we introduce the simplified notation 
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and insert our ansatz 03 into the DGLAP equation together with the LO expansion of 
the /9-function to get 
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Equating term by term in powers of L we find the recursion relation 


= (18) 
Po 

At this point we need to show that these recnrsion relations can be solved in terms of 
some initial condition and that they reproduce the exact LO solution in moment space. 
This can be done by taking Mellin moments of the recursion relations and solving the 
chain of these relations in terms of the initial condition Ao(x). At LO the solution of CO 
in moment space is simply given by 


= (-|-P(»>) q(N,a.(Q!)), (19) 

having imposed the initial condition Aq = g(x, aslQ^))- At this point we plug in this 
solntion into lO to obtain 
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which clearly coincides with m, after a simple expansion of the latter 
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Notice that this non-singlet solntion is an exact one. In the singlet case the same approach 
will sncceed at the same order and there is no need to introdnce trnncated solntion at this 
order. As expected, however, things will get more involved at higher orders, especially in 
the singlet case. 


The strategy that we follow in order to construct solutions of the DGLAP equations 
is all contained in this trivial example, and we can summarize our systematic search of 
logarithmic solutions at any order as follows: we 

1) define the logarithmic ansatz up to a certain perturbative order and we insert it 
into the DGLAP equation, appropriately expanded at that order; 

2) derive recurrency relations for the scale invariant coefficients of the expansion; 

3) take the Mellin moments of the recnrrency relations and solve them in terms of the 
moments of the initial conditions; 

4) we show, finally, that the solntion of the recnrsion relations, so obtained, is exactly 
the solntion of the original evolntion eqnation firstly given in moment space. This ap¬ 
proach is sufficient to solve all the equations at any desired order of accuracy in the strong 
coupling, as we are going to show in the following sections. Ultimately, the success of the 
logarithmic ansatz lays on the fact that the solntion of the DGLAP eqnations in QCD 
resnms only logarithms of the conpling constant. 
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4 Truncated solution at NLO. Non-singlet 


The extension of our procedure to NLO (non-singlet) is more involved, but also in this 
case proofs of consistency of the logarithmic ansatz can be formulated. However, before 
starting our technical analysis, we define the notion of “truncated solutions” of the DGLAP 
equations, expanding our preliminary discussion of the previous sections. We start with 
some definitions. 


A truncated solution retains only contributions up to a certain order in the expansion 
in the coupling. We could define a 1-st truncated solution, a 2-nd truncated solution 
and so on. The sequence of truncated solutions is expected to converge toward the exact 
solution of the DGLAP as the number of truncates increases. This can be done at any 
order in the expansion of the DGLAP kernels (NLO,NNLO,NNNLO,...). For instance, 
at NLO, we can build an exact solution in moment space (this is true only in the non¬ 
singlet case) but we can also build the sequence of truncated solutions. It is convenient to 
illustrate the kind of approximations which are involved in order to obtain these solutions 
and for this reason we try to detail the derivations. 


Let’s consider the NLO non-singlet DGLAP equation, written directly in moment 
space 
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Notice that equation is the exact NLO equation. In particular we have preserved 
the structure of the right-hand side, that involves both the beta function and the NLO 
kernels and is given as a ratio of two polynomials in Ug 
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is the NLO kernel. The factorization of the LO solution from the NLO equation can 
be obtained expanding the ratio P/jd in cks, which allows the factorization of a l/ag 
contribution. Equivalently, one can redefine the integral of the solution in moment space 
by subtraction of the LO part 


( P^^O{x,a) PLo{a) \ 
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Denoting by hi = /9i//9o, the truncated differential equation can be written as 
df{N,a,) 2 
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Notice that this solution of the truncated equation, exactly as in the exact solution (IP, 
contains as a factor the LO solution and therefore can be rewritten in the form 

(cKs — «o) f hi 


f{N, as) = exp 
where f^‘^{N,as) is given by 
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Eq. (1^ exemplifies a typical mathematical encounter in the search of solutions of PDE’s 
of a certain accuracy: if we allow a perturbative expansion of the defining equation 
arrested at a given order, the solution, however, is still affected by higher order terms in 
the expansion parameter (in our case a*). To identify the expansion which converges to 
jsni) proceeds as follows. We start from the 1-st truncated solution. 

Expanding JP to first order around the LO solution we obtain 
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/(TV, a,) = /^0(iV, a,) X 1 1 + } , (31) 

which is the expression of the 1-st truncated solution, accurate at order One can 
already see from ea that the ansatz which we are looking for should involve a double 
expansion in two values of the coupling constant: a* and cxq. This point will be made 
more clear below. Eor this reason we are naturally lead to study the logarithmic expansion 
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which is the obvious generalization of the analogous LO expansion US- 


Inserting this ansatz in the NLO DGLAP equation, we derive the following recursion 
relations for Ar, and 
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together with the initial condition 


f{x, QI) — ^0 + «s(Qo)-® 0 - 


(34) 


At this point we need to prove that the recnrsion relations dasi) reprodnce in moment 
space ea. To do so we rewrite the recnrsion relations in Mellin-space 

B„+i(iV) = -B„(N) - ^A„+,(]V) - |.p(»)(]V)B„(iV) - ^F(‘)(]V)>1„(]V), 

(36) 


and search for their solution over n. After solving these relations with respect to Aq and 
Bq, it is simple to realize that our ansatz (El exactly reproduces the truncated solution 
CD only if the condition i?o = 0 is satisfied. In fact, denoting by 
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the recursive coefficients, we can rewrite the recursion relations as 
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we identify the structure of the nth iterate in close form 
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Substituting the expressions for An and Bn so obtained in terms of Aq and Bq in the 
initial ansatz, and summing the logarithms (a procedure that we call “exponentiation “) 
we obtain 
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expression that can be rewritten as 
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This expression after a simple rearrangement becomes 
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This solntion in moment space exactly coincides with the trnncated solntion en if we 
impose the condition Bq = 0. It is clear that the solntion gets organized in the form of 
a donble expansion in the two variables ag and ckq. While ck^ appears explicitely in the 
ansatz (El, CKo appears only after the logarithmic summation and the factorization of the 
leading order solution. An obvious question to ask is how should we modify our ansatz 
if we want to reproduce the exact solution of the truncated DGLAP equation in moment 
space, given by eq. El. The answer comes from a simple extension of our recursive 
method. 


4.1 Higher order truncated solutions 


We start by expanding the solution of the truncated equation El> whose exponential 
factor is approximated by its double expansion in and ckq to second order, thereby 
identifying the approximate solution 
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To generate this solution with the recursive method it is sufficient to introduce the higher 
order (2nd order) ansatz 

f{x, «s) = X] — [An{x) + agBn{x) + alCn{x)] , (44) 

n=0 

where we have included some new coefficients Cn{x) that will take care of the higher order 
terms we aim to include. Inserted into the NLO DGLAP equation, this ansatz generates 
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an appropriate chain of recursion relations 
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that we solve by going to Mellin space and obtain 
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where the initial conditions are f(N, oq) = ^o and Bq = Cq = 0. It is a trivial exercise 
to show that the solution of the recursion relation, inserted into (in, coincides with (m> 
after exponentiation. 

Capturing more and more logs of the truncated logarithmic equation at this point is 
as easy as never before. We can consider, for instance, a higher order ansatz accurate to 
O(a^) for the NLO non-singlet solution 


^ jja 

f{x, aj = ^ — [An{x) + asBn{x) + alCn{x) + a^P„(a:)] , (47) 
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that generates four independent recursion relations. The relations for An+i, Pn+i, Cn+i 
are, for this extension, the same as in the previous case and are listed in (Uni). Hence, we 
are left with an additional relation for the D^+i coefficient which reads 
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These are solved in Mellin space with respect to Ao,Bq,Co, Dq (with the condition Bq = 
Cq = Dq = 0). We obtain 
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which is the solntion of the trnncated eqnation compnted with an 0 (q;^) accnracy. 


5 Non-singlet truncated solutions at NNLO 


The generalization of the method that takes to the trnncated solntions at NNLO is more 
involved, but to show the equivalence of these solutions to those in Mellin space one 
proceeds as for the lower orders. As we have already pointed out, one has first to expand 
the ratio P//5 at a certain order in a^, then solve the equation in moment space - solution 
that will bring in automatically higher powers of - and then reconstruct this solution 
via iterates. 


At NNLO the kernels are given by 

.2^NNL0 




and the equation in Mellin-space is given by 

df{N,as) 


da. 


pNNLO 
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We search for solutions of this equation of a given accuracy in ag and for this purpose we 
truncate the evolution integral of the ratio P/f3 to 0{al). This is the first order at which 
the component of the kernels appear. As we are going to see, this will generate the 
first truncate for the NNLO case. Therefore, while the first truncate at NLO is of 0{as), 
the first truncate at NNLO is of 0{al). We obtain 


I NNLO 


da 


’ Q^O 


Pnnlo{x, a) _ Plo{o) 
PNNLoiP) I3 lo{(R) 


-Pitto — —P2 Q(o a PiCKs + —R2a^g. (53) 
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At this retained accuracy of the evolution integral, the exact solution of the corre¬ 
sponding (truncated) DGLAP equation can be found, in moment space, as in (ITdl) 


f{N,a,) = f{N,ao) 


a. 


tto 


- 2 - 


=( 0 ) 


1 -|- (cTs — ao) 


p(l) ^ p(0)^^ 


7r/3o 27r/3^ 


+cr„ 


+“o 


p(l)2 p(2) p(0)p(l)^^ p{l)^^ P(0)^/?2 p{0)^2 p(0)^. 


2vr2/?2 

■p(l)2 


+ 


dvr^/do 27r2/3Q STr^/d, 

p(2) p(o)p(i)/?i p(i)/?i p(o)^/?2 


27 t ‘^( 3 q 47r^/9o 27r^/?, 


8vr2/?2 


87r2/?4 


167r2/3Q 

p(0)^2 

IbTT^/do 


+ 


167r2/3Q 

p^' 

Ibvr^/do 


+Q!ocrs 


p(i)2 


p{0)p(l)^^ p(0)2^2 


7r2/?3 




(54) 


At this order this solution coincides with the exact NNLO solution of the DGLAP equa¬ 
tion, obtained from an exact evaluation of the integral o, followed by a double epansion 
in the couplings. Therefore, similarly to (m, the solution is organized effectively as a 
double expansion in as and ao- This approach remains valid also in the singlet case, when 
the equations assume a matrix form. As we have already pointed out above, all the known 
solutions of the singlet equations in moment space are obtained after a truncation of the 
corresponding PDE, having retained a given accuracy of the ratio P/{3. For this reason, 
and to compare with the previous literature, it is convenient to rewrite dSlI) in a form 
that parallels the analogous singlet result [12]. It is not difficult to perform the match of 
our result with that previous one, which takes the form m, 121, Q 


f(N,a,) = t;(iV,a.)/io(iV,a.,a„)t/-'(]V,ao) 


H-oo 




K.= l 


fLo{N, Ug, ao) 


H-OO 


n -1 


1 + 5^ (7.(Af)aS 


K=1 


(56) 


which becomes, after some manipulations 


/(7V,a,)= ^ 

\ao 


, -J-p(O) 

as\ 1^0 


[l + {as-ao)U,{N) + ap2{N) 

-asaoUl{N) + al {Ul{N) - U^iN))] f{N, ao ), 


(56) 


where the functions Ui{N) are defined as 


U,{N) 

U2{N) 

R2{N) 

Ro{N) 


1 


6iPW(^) 


npo 

i [fif(]V) - B,(Af)] , 


pa)(iv) 


= Pi(iV), 


p(2)(JV) 

2vr2Ho 


^Pi(Ar) + 
dvr (dvr) 


;Ro{N) 




(57) 
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where /3i/^o = h, /? 2 /^o = h- 


We intend to show rigorously that this solution is generated by a simple logarithmic 
ansatz arrested at a specific order. For this purpose we simplify (jsni) obtaining 


f{N,Q^) = 


tto/ 


Ro{N) 


1 - i?i(iV)«o + ^Rl{N)al + ^R 2 {N)al + R,{N)as 


-Rl{N)asao + 


f{N,Ql), 


(58) 


and the ansatz that captures its logarithmic behaviour can be easily found and is given 
by 


fix, 0=)""“ 


E 

n=0 


An{x) 


n\ 


In 


OisiQ"^ 


O^siQl) 


n=0 




Cn(x) 


n=0 


n\ 


In 


+ o!s{q^) 

as{Q^] 
c^siQl 


BJx) 


n\ 


In 


C^s{Ql] 


(59) 


Setting the initial conditions as 


/(x, Ql) = Ao{x) + aoSo(x) + aoC'o(x), 


(60) 


and introducing the 3-loop expansion of the /3-function, we derive the following recursion 
relations 


An+l{x) = -^P®(x) ® An{x), 

Po 

Pn+l(x) = -Bn{x) - -^^An+l{x) - ^P^^\x) ® P„(x)- ]-P^^\x) 0 zl„(x), 

AttPo Po ttPo 

Cn+l{x) = -2Cn{x) - ^^Bn{x) - ^Pn+l(x) - ^^^/l„+i(x) 

-^P^°\x) ® Cn(x) - ]-P^^\x) ® Bn{x) 

Po 'xPo 


We need to show that the solution of the NNLO recursion relations reproduces JSEl) in 
Mellin-space, once we have chosen appropriate initial conditions for 74o(W), Bq{N) and 
C,{N) \ _ 

^It can be shown that the infinite set of recursion relations have internal symmetries and different 
choices of initial conditions can bring to the same solution. The choice that we make in our analysis is 
the simplest one. 
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At NLO we have already seen that Bo{N) has to vanish for any N, i.e. Bq{x) = 0, 
and we try to impose the same condition on Co{N). In this case we obtain the recnrsion 
relation for the moments 

C„MN) = -2C„(N) - 

_|_p{0)(iv)C'„(iV) - ^pd)(iv)P4iV) 

Po 7r/to 

-^P'^\N)A„(N), ( 62 ) 

which combined with the relations for An{N) and Bn{N) give 

= I -Rl (i?o - 1)" - Ir2R^ + ^ [(i?o - 2)" Rj + {Ro - 2)" i?2 + RlR^] } 

Bn = [K - (^0 - 1)”] RlAo, 

An = i?o Ao , (63) 

where the N dependence in the coefficients Ri has been snppressed for simplicity. The 
solntion is determined exactly as in (EH) and (EHI) and can be easily bronght to the form 

1 /a 

f{N,Q'^) = - ( — j ^2 — 2Ri{as — ao) + Rl{as — O'oY A R2{oil — Oil)] Aq, (64) 

2 V«o/ 

which is the result quoted in eq. JSHI) with Ao = f{N,ao). We have therefore identified 
the correct logarithmic expansion at NNLO that solves in x-space the DGLAP equation 
with an accuracy of order 


6 Generalizations to all orders : exact solntions of trnn- 
cated eqnations bnilt recnrsively 

We have seen that the order of approximation in of the truncated solutions is in 
direct correspondence with the order of the approximation used in the computation of the 
integral on the right-hand-side of the evolution equation jP. This issues is particularly 
important in the singlet case, as we are going to investigate next, since any singlet solution 
involves a truncation. We have also seen that all the known solutions obtained in moment 
space can be easily reobtained from a logarithmic ansatz and therefore there is complete 
equivalence between the two approaches. We will also seen that the structure of the ansatz 
is insensitive to whether the equations that we intend to solve are of matrix forms or are 
scalar equations, since the ansatz and the recursion relations are linear in the unknown 
matrix coefficients (for the singlet) and generate in both cases the same recursion relations. 

We pause here and try to describe the patterns that we have investigated in some 
generality. We work at a generic order N^LO, with m=l denoting the NLO, m = 2 the 
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NNLO and so on. We have already seen that one can factorize the LO solution, having 
defined the evolution integral 


j I ^ ^ fPN^Lo{x,Oi) Plo{<^) 

^N^Lo{Ois,ao) — / da \ — -——-———— 

V PN^Loia) PLo[a) 


' OiQ 


and the exact solution can be formally written as 

f{N,as) = f°{N,as) x 


(65) 


( 66 ) 


A Taylor expansion of the integrand in the jnsi) around = 0 at order k = {m — 1) gives, 
in moment space. 


/ Qg) 

V ^N^Loiois) 


^ Ri{P^^\P^^\ N) + R2{P^^\ N)as 

PLOyO^s) J 

+R^{P^^\P^^\P^^\N)al 

+ ... + /?,+i(p(°), P('),..., P("^), iVX, (67) 


which at NNLO becomes 


/ Pnnlo{N, a^) 

V ^NNloiois) 


PLo(as) A 

Phoiois) J 


^ Ri{P^°\P^^\N) + R2{P^^\P^^\P^^\N)as.{68) 


Thus, integrating between as and oq we, obtain the following expression for lN-^Lo{c>is, ao) 
at 0{ag) 

^N^Loi^s, Cto) = —Pl«0 ~ o-^2«0 ~ • • •-Pk«o + Rl^^s + -P2«s + • • • H- RnOig , (69) 

Ik Ik 

where the (P^°\ pP\ ..., p(”^), N) dependence in the P^ coefficients has been omitted. 
To summarize: in order to solve the K-truncated version of the N'^LO DGLAP equation 


df{N, as 
das 


Pn">lo(A^, ds) 

/^N“Lo(tts) 


/(iV,«.), 


(70) 


which is obtained by a Taylor expansion - around a* = 0 - of the ratio Pn™lo(-^) «s)//^N™Lo(cts)) 
we need to solve the equation 


df{N, as) 

das 


— [Po + a^Pi + Q:gP2 + ... + a^PA f{N,as), 

r\i -* 


(71) 


where the coefficients P^ have a dependence on P^^'> and P^d in the NLO case, and on 
P(°), Pd) and P(2) in the NNLO case. Eq. ea admits an exact solution of the form 


f{N,as)={^ 
ao 


Rq 


exp 


Pi («s ~ «o) + X-R2 (q^s ~ '^0) + • • • H- Rk (c^s ~ '^0) f /(-^I “o)) 

Z K 1 


( 72 ) 


having factorized the LO solution. 
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At this stage, the Taylor expansion of the exponential around (Q!s,tto) = (0; 0) gener¬ 
ates an expanded solution of the form 


(n) 


fNmLo{N,as) ^ fLoiN,ao)e^^^^o 


- fioiN, ao) ( 1 + In^lo + ^ (^^N^Lo) + • • • 


(73) 


which can be also written as 


fLo{N, ao) + lifliLo + 2! + • • • ) — 

CTs) X [co + CHg (ci^o + Ci^iCTo + Ci,2tto + • • • + Ci^k-I^o ^ + Ci^^O^q) 

+<^5 (c2,0 + C2 ,iQ!o + 02,2^0 + • • • + C2,Ktto) + • • • -!-■■■], 

(74) 

where the coefficients Cij are defined in moment space. The K-th truncated solution of the 
equation above 


/w-Lo(^,a.)lo(aj) = f^°iN,as) I 

\i,j=0 




CX^OiQCij I , 


(75) 


is therefore accurate to O(a^), and clearly does not retain all the powers of the coupling 
constant which are, instead, part of JZHl). However, as far as we are interested in an 
accurate solution of order k, we can reobtain exactly the same expression from a:-space 
using the ansatz 


fw^Loix, as)lo(aj') 


(^^n(2^) + -f a‘^A‘^(x) + . . . -1- CKg (x)^ 

n=0 



(76) 

which can be correctly defined to be a truncated solution of order k! of the («:)-truncated 
equation. Since we have truncated the evolution integral at order k, this is also the 
maximum order at which the truncated logarithmic expansion JZEl) coincides with the 
exact solution of the full equation. This corresponds to the choice n' = k. Notice, however, 
that the number of coefficients A'^ that one introduces in the ansatz is unrelated to k and 




can be larger than this specific value. This implies that we obtain an improved accuracy 
as we let k' in dZEI) grow. 


If we choose the accuracy of the evolution integral to be k, while sending the index 
k' in the logarithmic expansion of JZED to infinity, then the ansatz that accompanies this 
choice becomes 


OO 

fN^Lo{N,as) = 'Y 

n=Q 




. 1=0 



In 


WiQDJ 


n 




(77) 
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and converges to the exact solution of the (order k) truncated equation (TTHll . Obviously, 
this exact solution starts to differ from the exact solution of the exact DGLAP equation 
at Also in this case, as before, one should notice that the double expansion 

in Us and ao of the exact solution of the (K)-truncated equation can be reobtained after 
using our exponentiation, and not before. We remark, if not obvious, that the recursion 
relations, in this case, need to be solved at the chosen order k!, as widely shown in the 
examples discussed before. 

We remark that, as done for the LO, we could also factorize the NLO solution and 
determine the N'^LO solution using the integral 


In^lo 



f Pn”^lo{x, a) 

V (^N-^Lo{0i) 


Pnlo{c^) a 

Pnlo{(P) J ’ 


(78) 


and then restart the previous procedure. Obviously, the two approaches imply a resum¬ 
mation of the logarithmic behaviour of the pdf’s in the two cases. 

It is convenient to summarize what we have achieved up to now. A truncation of 
the evolution integral introduces an approximation in the search for solutions, which is 
controlled by the accuracy (k) in the expansion of the same integral. The exact solution 
of the corresponding truncated equation, as we have seen from the previous examples, 
involves all the powers of ag and ao and, obviously, a further expansion around the point 
tts = ao = 0 is needed in order to identify a set of truncated solutions which can be 
reobtained by a logarithmic ansatz. This is possible because of the property of analiticity 
of the solution. Therefore two types of truncations are involved in the approximation of 
the solutions: 1) truncation of the equation and 2) truncation (k') of the corresponding 
solution. In the non-singlet case, which is particularly simple, one can therefore identify 
a wide choice of solutions (by varying n and k') that retain higher order effects in quite 
different fashions. Previous studies of the evolution using an ansatz a la Rossi-Storrow 
pm, borrowed from the pdf’s of the photon, were therefore quite limited in accuracy. 
Our generalized procedure is the logical step forward in order to equate the accuracy of 
solutions obtained in moment space to those in x-space, without having to rely on purely 
numerical brute force methods. 

3 


6.1 Recursion relations beyond NNLO and for all /t’s 

In the actual numerical implementations, if we intend to use a generic truncate of the 
non-singlet equation (the result is actually true also for the singlet), it is convenient to 

^In (21 a flag variable called IMODEV allows to switch among the exact solution (IMODEV=l), the 
exact solution of the 0(ag) truncated equation (IMODEV=2). A third option (IMODEV=0) involves at 
NLO and at NNLO 0(as) and Ofog), respectively, truncated ansatz that, in our cases, are reconstructed 
logarithmically. 
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work with implementations of the recursion relations that are valid at any order. In fact 
it is not practical to rederive them at each new order k' of approximation. Next we are 
going to show how to do it, identifying generic relations that are easier to implement 
numerically. 


The expression to all orders of the DGLAP kernels is given by 

z+i 




da. 


(91ng2 


1=0 

OO 

■E 

i=0 


a 


i+2 


A 


[An] 


i+i • 


and the equation in Mellin-space in Melin space is given by 

\ / + 1 


dAN.a.) _ TlZaiHy P^'HN) 


da. 


r.i+2 Ai 

5 ( 477 )^+^ 


a* 

Z^i=0 


f(N,a,), 


(79) 


(80) 


whose exact solution can be formally written as 


. _^p(0) 

«o/ 


(81) 


where JF is obtained from the evolution integral and whose specific form is not relevant 
at this point. 

We will use the notation (P, P) to indicate the sequence of components of the kernels 
and the coefficients of the /9-function (P®, p9), p(2) ) 

Then, the Taylor expansion around o;* = ckq of the solution is formally given by 


^T{as,ao,P,0) _ 


J2^n{dJ^,d^J^, ...,d^J^)U=^Pas - aoT^^ 


(82) 


n=0 


for an appropriate <hn(i9P, (9^P,..., d'^yp). <!>„ is a function that depends over all the partial 
derivative obtained by the Taylor expansion. Since it is calculated for o;* = ckq, it has a 
parametric dependence only on ao, and we can perform a further expansion around the 
value ckq = 0 obtaining 


..fKpoPA ^ 0^0 ^ 

^ m\ daP" 

m=0 ^ 


<hn(ao, P, P) {ois - ao)” -r 

n\ 


n=0 


oo=0 


(83) 


This expression can always be arranged and simplified as follows 


^T{as,ao,P,l3) _ 


E 

n=0 


a' 


n\ 


a 


,n —1 


a' 


P, /9) + -n^$A0, P, P) + A P) 


n\ 


n\ 


«o=0 


ao+ 


n—2 ^n—1 

$^0, P, P) - n^d<l>n{ao, P, P) 

n\ 


1 (n — l)nQ;'; " ^ a 


nl 


oo=0 


20 



A/3) 


«o=0 


“o + (• • O^o + 


l + a,d°’(A« + ... + <?f(A/3)) 


+ao (^^‘’(^J) + a4'\pj) + aUi'-\PJ) + ... + 

+ol (i'i\PJ) + a£\pj) + all^i\pj) + ... + a’Xf>{Pj)) 


+a« {(ir\P, P) + P) + aUt’(Pj) + ■ ■ ■ + a:(i"'HP,P)) 

oo oo 

n=0 m=0 

(84) 

where we are formally absorbing all the dependence on both the kernels P and the /3- 
function f3, coming from the functions 5™$^ calculated at the point ao = 0, in the 
coefficients Finally, we can reorganize the solution to all orders as 

^ p(0) oo OO 

-) ” /(iV.ao)EE“S‘“"««’"’(^"'3)- (85) 

^ ^ n=0 m=0 

With the help of the general notation 

Pnnlo = {P^''\P^^\P^"^), 

Pnlo = (/5o,/3l), 

Pnnlo = {Po,Pi,P2 ), (86) 

we can try to indentify, by a formal reasoning, the exact solution up to a hxed - but 
generic - perturbative order of the expansion of the kernels. 

To obtain the NLO/NNLO exact solution it is sufficient to take as null the components 
(p(2)^ p(3)^ (^^ 2 , /? 3 ,...) for the NLO and {P^^\ ...) and (/Sa,...) for the NNLO case. 

Since eq. JHHI) contains all the powers of aoUs up to (i.e. it is a polynomial expression 

of order if we aim at an accuracy of order a^, we have to arrange the exact 

expanded solution as 

/(Va.)=(- ° m,a„)Y.Y,aia:-ii<P(Pj)+0(a^), (87) 

with 0 (q(^) indicating all the higher-order terms containing powers of the type ckoCk^ -1- 
... -|- CDoCK^. Hence, the index k represents the order at which we truncate the solution. 
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As a natural generalization of the cases discussed in the previous sections, we introduce 
the higher-order ansatz (/t-truncated solution) 


m a.) = Y, 

n=0 


m=0 


0”(iV) 


-a. 


n\ 



( 88 ) 


that reproduces the exact solution JHIl) expanded at order k. 

In fact, inserting this last ansatz in (IHID we generate a generic chain of recursion 
relations of the form 


0;+i(/V) = F»(0:(pl"l,/3„)), 

0;+i(iV) = F‘(0.”,0^i,Oi,Pj), 

= F”(0»,...,0.«,0^i,...,0-!,F,,3) . (89) 

A deeper look at the explicit structure of JEHI), for the NLO non-singlet case, reveals 
the following structures for the generic iterates 


0."+i(A') = -|;[i’'"’(A')0:(/V)], 

0;:+i(A') = ~ [F'">03 (N) - ^ [P"-\N)Ol-\N)] 

- J^O-i(Af) - nO:(N) -(k- 1) J^Or'(A'), (90) 

at NLO, while for the NNLO case we obtain 


0°+i(iV) 

OlAN) 


= --[pW(iV)0°(iV)], 

= [pW(iV)Oi(iV)] - ^ [pW(iV)O0(iV)] 


A 


47r/3o 


Ol^,{N)-Ol{N), 


P [Fl»)(Af)0;(]V)] - ^ [F<‘)(Ar)Or‘(iV)] 


27r2/?o 


[p(2)(Ar) o;;-2(iV)] 


A 


47r/3o 


QK-l 


(iV)- 


A 


IGtt^/So 


qk-2 

'^n+1 


(AT) 


A 


kO:{N) -{k- l)^or^(iv) -{k-2) 


A 




or^(iv). 


(91) 
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Hence, one is able to determine the structure of the K-th recursion relation when the k = 0 
and K = 1 cases are known. This property is very useful from the computational point of 
view. ^ 

Passing to the resolution of the recursion relations in moment space, we get the formal 
expansion of each O^(iV) in terms of the initial condition Oo(A^), which reads 

/ {N, ao) = Oq{N) + 0 tto + 0 ttg + ••• + 0 «o • (92) 

Here we have set to zero all the higher order terms, as a natural generalization of Bq = 
Co = 0 ..., according to what has been discussed above. 

These relations can be solved as we have shown in previous examples, and the generic 
structure of their solution can be identified. If we define Rq = — the expressions 

of all the 0!^{N) in terms of f{N, ao) = /o become 

0°(iV) = P-/o 

Ol{N) = G\Rl{Ro-ir,Pj)fo, 

Ol{N) = G2(Po", (^0 - 1)", (i?o - 2)", P, /3)/o, 


Ol{N) = G^{Rl (Po - 1)”, {Ro - 2)",..., (Po - P, P)fo , 


(93) 


and in particular, by an explicit calculation of 0'^{N), one can work out the structure of 
these special functions G^. For instance, we get for m = k the expression 

K, 

G'{Tr„ (R, - 1)”, (R - 2)“,.... (Rj - k)”, = ir&iP, n . (94) 

j=0 

for suitable coefficients Substituting the 0^{N) functions with m = 0,... k in the 
higher-order ansatz (IHHD and performing our exponentiation, we get an expression of the 
form 



(95) 


which can be written as follows by the use of eq. (ini 

( \ Ro Km 

- (96) 

^ m=0 j=0 

This is the exact solution expanded up to 0{ag) in accuracy. 

^The relations pi|l and m hold also in the NLO/NNLO singlet case and can be generalized to any 
perturbative order in the expansion of the kernels. 
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7 The Search for the exact non-singlet NLO solntion 


We have shown in the previous sections how to construct exact solutions of truncated 
equations using logarithmic expansions. We have also shown the equivalence of these ap¬ 
proaches with the Mellin method, since the recursion relations for the unknown coefficient 
functions of the expansions can be solved to all orders and so reproduce the solution in 
Mellin space of the truncated equation. The question that we want to address in this 
section is whether we can search for exact solutions of the exact (untruncated) equations 
as well. These solutions are known exactly in the non singlet case up to NLO. It is not 
difficult also to obtain the exact NNLO solution in Mellin space, and we will reconstruct 
the same solutions using modified recursion relations. The expansions that we will be 
using at NLO are logarithmic and solve the untruncated equation. The NNLO case, in¬ 
stead, will be treated in a following section, where, again, we will use recursion relations 
to build the exact solution but with a non-logarithmic ansatz 

The exact NLO non-singlet solution has been given in eq. (121. The identification of 
an expansion that allows to reconstruct in moment space eq. (121 follows quite naturally 
once the typical properties of the convolution product ® are identified. For this purpose 
we define the serie of convolution products 

n\ 

n=0 

that acts on a given initial function as 

- x) ®+FaPa{x) +^FIPa® Pa<^ +■ 

= (j){x) + Fa {Pa ® 0) + ^F^ {Pa ® Pa<^) 0(a:) -k .... (98) 

The functions Fa and Fb are parametrically dependent on any other variable except 
the variable x. The proof of the associativity, distributivity and commutativity of the ® 
product is easily obtained after mapping these products in Mellin space. For instance, for 
generic functions a{x) b{x) and c{x) for which the ( 8 ) product is a regular function one has 

M [(a ® 6 ) 8 ) c] {N) = M[a®{h® c)] {N) 

= a{N)B{N)C{N), (99) 

where A4 denotes the Mellin transform and N is the moment variable. Also one obtains 

^FaPa{x)® gPsPflW® 0 ( 3 .) = ^{FaPa{x)+FsPb{x))® ( 100 ) 

and 

M 0] {N) = (101) 

®In PEGASUS |2l the NNLO non-singlet solution is built by truncation in Mellin space of the evolution 
equation, while the NLO solution is implemented as an exact solution. 
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with M ^-independent, since both left-hand-side and right-hand-side of (USD can be 
mapped to the same fnnction in Mellin space. Notice that the role of the identity in 
®-space is taken by the fnnction 5(1 — x). We will also nse the notation 

(j){x) = {Aq{x) +FaAi{x) ®-\-.. .)(j){x), ( 102 ) 

\n=0 / g, 

where the A'n{x) and the An{x) capture the operatorial and the functional expansion - 
respectively - and are trivially related 

A'^(x) ® 4>{x) = An{x). (103) 

To identify the x-space ansatz we rewrite (1231) as 


/(TV, a) 




/(A^,«o) 





where we have introduced the notations 


(104) 


L 

M 

a{N) 

b{N) 


1 

log—, 

ao 


= log 


+ OisPi 


4:71 Pq -|- OqPi ’ 
2P^^\N) 


/3o ’ 

2p(o)(iV) 4P(i)(iV) 


/^o 


Pi 


Our ansatz for the exact solution in x-space is chosen of the form 




vn=0 
oo s 


n\ 




\rn=0 


ml 


f{x,Ql 




5=0 n =0 

CO s 

TT 


n!(s — n)\ 


Q(^) 


n i\/rs-n 


L^M 


(105) 


( 106 ) 


where in the first step we have turned the product of two series into a single series of a 
combined exponent s = n + m, and in the last step we have introduced the functions 


Cn{.x) = A'^ix) ® P'_„(a;) ® /(x, Ql), (n < s). 


(107) 
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Setting Q = Qo in (fTTlHll we get the initial condition A'q{x) = Bq{x) = 5(1 — x) or, 
eqnivalently, 

f{x,Ql) = C^{x). (108) 

Inserting the ansatz (Tmoli into the NLO DGLAP eqnation, with the expressions of the 
kernel and beta function included at the corresponding order, we obtain the identity 


EE 



Pi 2 
a 


Ibvr^ 


r^s+i 

'^n+l 


Pi 2^, 


s+1 


27r ^ 4j^2 


a 




n!(s — n)\ 
n!(s — n)\ 


(109) 


Equating term by term the coefficients of a and we find from this identity the new 
exact recursion relations 


Qt! = (110) 

cy' = -C‘„%\-JPP>(SC‘„, (111) 

or, equivalently, 

C‘^ = (112) 

C‘n = -CU, - jP"^ K C‘-K (113) 

Notice that although the coefficients G® are convolution products of two functions, the 
recursion relations do not let these products appear explicitely. These relations just 
written down allow to compute all the coefficients G® {n < s) up to a chosen s starting 
from Gq, which is given by the initial conditions. In particular eq. (HH allows us to 
move along the diagonal arrow according to the diagram reproduced in Table HJ eq. (I1 1 311 
instead allows us to compute a coefficient in the table once we know the coefficients at 
its right and the coefficient above it (horizontal and vertical arrows). To compute G® 
there is a certain freedom, as illustrated in the diagram. For instance, to determine G| 
we can only use CH, and for the coefficients Gq we can only use nn. For all the other 
coefficients one can prove that using dm or nn brings to the same determination of 
the coefficients, and in our numerical studies we have chosen to implement dm, being 
this relation less time consuming since it involves P® instead of 

The recursion relations defining the iterated solution can be solved as follows. From 
the first relation dini), keeping the s-index fixed, we have 


G® = 


_Ap(o)^ 

Po 

Po 


/OS —1 

^n-1 


C s—n—1 
0 


(114) 
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then, since the second relation (llld^l also holds for n = 0, we can write (using eq. (11141) 1 




4 

Ji 


Co = 


C^o = 


^piO) _ Ap(1) 

iPo A 

^.p(0) _ Ap(l) 


C S—1 
0 


_Po Pi 

Finally, inserting the above relation in (HH we can write 


Oq. 




|.p(01 

n 

|.p(0) _ 4 pdl 

_Po 


_Po Pi 


/^O 
'-'0 ’ 


( 116 ) 


(116) 


which is the solution we have been searching for. The last step in the proof consists in 
taking the Mellin transform of this operatorial solution and summing the corresponding 
series 


f(N,a,) = XIE 


c;(iv) 


n nyTS-n 


L^M' 


= EE 

s=0 n=0 

that after summation gives 
/(/V,a.)=e-i^‘“’<"'‘“'(5S) 




o' 

CM I' 
1_ 

n 

|.p(0) _ i.p(i) 

_Po 


_Po Pi 


C^(N), (117) 


exp 


^pm(jv)-lpW(N) 




( 118 ) 


which is exactly the expression in eq. uni. Hence, it is obvious that the exact solution 
of the DGLAP equation (1221) can be written in x-space as 




(119) 


therefore proving that the ansatz (fTTibll reproduces the exact solution of the NLO DGLAP 
equation from x-space. 

A second version of the same ansatz for the NLO exact solution can be built using 
a factorization of the NLO DGLAP equation. This strategy is analogous to the method 
of factorization for ordinary PDF’s. For this purpose we define a modified LO DGLAP 
equation, involving 


dfLo{x,as) 

das 


Ols \ 
2 ttP^^^ ) 


P®(x) ® fLo{x,as), 


( 120 ) 
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whose solution is given by 


2p(^) 

/^O 


fLo{x,a) = e 
fLo{x,a) = e 


M 


/30 


fLo{x,a) 

S) 

f(x,ao), 


( 121 ) 


where we have introduced the ordinary LO solution fio, expressed in terms of a typical 
initial condition f{x,ao), and the NLO recursion relations can be obtained from the 
expansion 


fNLo{x,a) = 


BJx) 


vn=0 


n\ 


fLo{x,a). 


Inserting this relation into we obtain the recursion relations 


Bo{x) = (5(1-x), 
which is solved in moment space by 


B„ 




Bo{N) 


Bn{N) = 

Bo{N) = 1. 

The solution eq. (II 2211 can be re-expressed in the form 


( 122 ) 


(123) 


(124) 


fNLo{x, a) 



fix,ao) 


M 

e 



_2p(0) 5 
/30 J 


f{x,ao) 


(125) 


which agrees with ()l()4jl once a{N) and b{N) have been defined as in (|1()5|) . 

We have therefore proved that the exact NLO solution of the DGLAP equation can 
be described by an exact ansatz. Since the ansatz is built by inspection, it is obvious that 
one needs to know the solntion in moment space in order to reconstrnct the coefficients. 
Thongh the recnrsive scheme nsed to constrnct the solntion in x-space is more complex 
compared to the recursion relations for the truncated solution, its numerical implemen¬ 
tation is still very stable and very precise, reaching the same level of accuracy of the 
traditional methods based on the inversion of the Mellin moments. 


8 Finding the exact non-singlet NNLO solution 

To identify the NNLO exact solntion we proceed similarly to the NLO case and start from 
the DGLAP equation in moment space at the corresponding pertnrbative order 
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i \ 

a - 

i \ \ 

/^2 ^ ^2 
*-^0 '^1 *-^2 

i \ \ \ 

/^3 ^_ /^3 /^3 /^3 

*^0 '-^1 '^2 '-^3 

i \ \ \ \ 


Table 1: Schematic representation of the procedure followed to compute each coefficient 

Q- 


a/(iV,a,) _ (^)p(o)(iV) + (f)^P«(iV) + (f)=^p(2)(iV) 


da^ 




Pi 


■a^ + 


(^2 


az 


f{N,a.). 


(126) 


After a separation of variables, all the new logarithmic/non logarithmic and dependences 
come from the integral 

(127) 


'o‘s{Qo) 


/3 


NNLOi 


a 


and the solution of (innD is 


f{N,a) = /(AT, 


-M^ 


vn=0 


\m=Q 


ml 
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c{Ny 
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QP , (128) 


where we have defined 


C 

M 

Q 

a{N) 

b{N) 

c{N) 


167r^/?0 + dTTCK/di + Q!^/d2 

IfivT^/^o + dvrcKo/di + alP2 ’ 

1 27r(Q( - gp) 4/3 o/? 2 - /?i 

\/4:PqP2 — 27r(87r/3o + (a + «o)/3i) + Q(Q(o/32 

2PW(Ar) 

p(o)(jv) 4P(2)(Ar) 

Po P2 

‘^pWf^N) - 8P^^\N) + ^p(2)(Ar). 

Po P2 


(129) 

(130) 

(131) 

(132) 

(133) 

(134) 


Notice that for nf = 6 the solution has a branch point since 4/3o/32 —/^i < 0. If we increase 
Uf as we step up in the factorization scale then, for nj = 6, Q is replaced by its analytic 


29 


















continuation 


Q 


— , =arctanh 

- W 2 


27r(a - gp) \/- 4/jo^2 
27r(87r/3o + (a + ao)/3i) + aao/32' 


(135) 


Eq. (11281) incorporates all the nontrivial dependence on the coupling constant a* (now 
determined at 3-loop level) into £, M and Q. 


As a side remark we emphasize that it is also possible to obtain various NNLO exact 
recursion relations using the formalism of the convolution series introduced above. For 
this purpose it is convenient to define suitable operatorial expressions, for instance 


^2 



■JNNLUJ^ 


•> 





(136) 


which are manipulated under the prescription that the integral in a is evaluated before 
that any convolution product acts on the initial conditions. The re-arrangement of these 
operatorial expressions is therefore quite simple and one can use simple identities such as 
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Jnnlo 
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27ry 
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2 p(i)(x, a) 

27rz 
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fjNNLO^^^ 


Q 2 


^p(o) (g) 
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Po ^ Po ^ Po ^ ^ 

-- -8QP^^%, 


Po Po 


(137) 


to build the NNLO exact solution using a suitable recursive algorithm. For instace, using 
unni) one can build an intermediate solution of the equation 


given by 


dfNLo{x,as) P^^^ 


da. 


pNNLO 


fNLo{x, as 


(138) 


JnLO — 


(139) 


and then constructs with a second recursion the exact solution 

f{x,as) = S 2 fNLO- (140) 

A straightforward approach, however, remains the one described in the previous section, 
that we are going now to extend to NNLO. In this case, in the choice of the recursion 
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relations, one is bound to equate 3 independent logarithmic powers of £, A4 and Q that 
appear in the symmetric ansazt 


/(X,C3^) = 


and where 




n\ 

\n=0 

oo s t 


\m=0 


ml 


fMl) 


\P=0 


p\ 




5=0 t =0 n =0 
oo s t 

SSS 


n\{t — n)\{s — t)\ 
Dln{x) 




Dlnix) = A'^{x) ® (g) C'_i(a;) g) /(x, Ql). 


(141) 


(142) 


The ansatz is clearly identified quite simply by inspection, once the structure of the 
solution in moment space (II2811 is known explicitely. In (114111 we have at a first step 
re-arranged the product of the three series into a single series with a given total exponent 
s = n + m + p, and we have introduced an index t = n + m. The triple-indexed function 
Dl^[x) can be defined also as an ordinary product 


A%(^) = Anix) g Bt-n{x) g Cs-t{x), {n<t <s), 
where we have absorbed the g operator into the definition of A, B and C, 
A{x)B{x)C{x) = A'{x) g (^B'{x) g (^C'{x) g f{x, Ql) 
Setting Q = Qo in (HUD we get the initial condition 

f(x,Ql) = Dl„(x). 


(143) 


(144) 


(145) 


Inserting the ansatz (HUD into the 3-loop DGLAP equation together with the beta func¬ 
tion determined at the same order and equating the coefficients of a, and a^, we find 
the recursion relations satisfied by the unknown coefficients D^^{x) 


or equivalently 




jjs+1 
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s+i 

t,n 


= -2ft (A'S„ + A‘i‘„+i) - 8P“> ® A'„. 


A‘„ = -^P'“’®Aftgi. 


= 




Din = -2/3i (A+l„ + P,■+!,„+.) - 8P''> ® A.y ■ 


(146) 

(147) 

(148) 

(149) 

(150) 
(161) 
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n4 
-^4,4 

n4 n4 

^3 3 -^4 3 

n4 ^4 n4 

■^2 2 "^3 2 "^4 2 

^ ^ ^ ^ 

/ / / i 

n4 ^_ n4 ^_ n4 ^ 

4^1,0 ^ 4^2,0 ^ 4^3,0 ^ 4^4,0 

Table 2: Schematic representation of the procednre followed to compnte each coefficient 
Dl^ for s = 4. The nnderlined coefficients are compnted via eq. (ITIqI) . 

In the compntation of a given coefficient -Df^, if more than one recnrsion relation is 
allowed to determine that specihc coefficient, we will choose to implement the less time 
consuming path, i.e. in the order cini), (nnn) and (nnni). At a hxed integer s we proceed 
as follows: we 

1 . compute all the coefficients Dl^ with n 7 ^ 0 using (fTlQll : 

2 . compute the coefficient D^q using dnni); 

3. compute the coefficient with t ^ s using jnn), in decreasing order in t. 

This computational strategy is exemplihed in the diagram in Tabled for s = 4, where the 
various paths are highlighted. 

Following a procedure similar to the one used for the NLO case, we can solve the 
recursion relations for the NNLO ansatz with the initial conditions Zi)Qo(^)- Solving the 
relations (11491113111 . we obtain the chain conditions 
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Then, from the second relation we get the additional ones 
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(162) 


(153) 
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From the last relation we also obtain the relations 


A'„ = - 2 A - SP‘" 
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(154) 


which solve the recnrsion relations in x-space in terms of the initial condition Hq q. Finally, 
the explicit expression of the „ coefficient will be given by 
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(155) 


The solntion of the NNLO DGLAP eqnation reprodnced by Hill) in Mellin space will 
then be written in the form 
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( 166 ) 


which is equivalent to 

/()V,a,) = J-i^‘“’]'°dSj)exp 
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DUN), (157) 


_ p2 27r(87r/lo + {as + ao)Pi) + asaoP2, 

and it reproduces the result in (inHi). In x-space the above solution can be simply written 


as 


f{x,as{QU = e 




^ exp < log 
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(158) 


We have shown how to obtain exact NNLO solutions of the non-singlet equations using 
recursion relations, ft is clear that the solution shown above conceals all the logarithms 
of the coupling constant into more complicated functions of Ug and therefore it performs 
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an intrinsic resummation of all these contributions, as obvious, being the exact solution of 
the non-singlet equation at NNLO. A numerical implementation of the recursion relations 
associated to these new functions of the coupling constants, in this case, is no different 
from the previous cases, when only functions of the form log(Q;/Q;o) have been considered, 
but with a faster convergence rate. 


9 Truncated solutions at LO and NNLO in the singlet 
case 

The proof of the existence of a valid logarithmic ansatz that reproduces the truncated 
solution of the singlet DGLAP equation at NLO is far more involved compared to the non 
singlet case. Before we proceed with this discussion, it is important to clarify some points 
regarding some known results concerning these equations in moment space. First of all, 
as we have widely remarked before, there are no exact solutions of the singlet equations in 
moment space beyond those known at LO, due to the matrix structure of the equations. 
Therefore, it is no surprise that there is no logarithmic ansatz that can’t do better than 
to reproduce the truncated solution, since only these ones are available analytically in 
moment space. If we knew the structure of the exact solution in moment space we could 
construct an ansatz that would generate by recursion relations all the moments of that 
solution, following the same strategy outlined for the non-singlet equation. Therefore, 
inverting numerically the equations for the moments has no advantage whatsoever com¬ 
pared to the numerical implementation of the logarithmic series using the algorithm that 
we have developed here. However, we can arbitrarily improve the logarithmic series in 
order to capture higher order contributions in the truncated solution, a feature that can 
be very appealing for phenomenological purposes. 

The proof that a suitable logarithmic ansatz reproduces the truncated solution of the 
moments of the singlet pdf’s at NLO goes as follows. 

9.1 The exact solution at LO 

We start from the singlet matrix equation 

d f q^^\x,Q^) \ Pqq{x,as{Q^)) Pgg{x, ^ q^P{x,Q^) \ 

(91ogQ2 g{x,Q^) ) \ Pgq{x,as{Q^)) Pgg{x,as{Q^)) ) \ gix^Q"^) 

(159) 

whose well known LO solution in Mellin space can be easily identified 

/(iV, as) = L{as, tto, N)f{N, ao), (160) 

/ xRo(JV) 

and where L(as, ag, -^) = ( ^ ) is the evolution operator. 
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Diagonalizing the Rq operator, in the equation above, we can write the evolution 
operator ao, N) as 

r+ 


L{as,ao,N) = e+ 


a. 


ao 


+ e_ 




ao 


( 161 ) 


where r± are the eigenvalues of the matrix Rq = {—2/Po)Po and e+ and e_ are projectors 
ns, nsi defined as 

1 


e± = 


Ro - r^I 


r± — Tzsp L 

Since the e± are projection operators, the following properties hold 


( 162 ) 


e+e+ = e+ e_e_ = e_ 

Hence it is not difficult to see that 


e+e_ = e_e+ = 0 


Ro{N) = e+r+ + e_r_ . 


e+ + e_ = 1. (163) 


(164) 


It is important to note that one can write a solution of the singlet DGLAP equation in a 
closed exponential form only at LO. 

It is quite straightforward to reproduce this exact matrix solution at LO using a 
logarithmic expansion and the associated recursion relations. These are obtained from 
the ansatz (here written directly in moment space) 


subject to the initial condition 
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n=0 


n\ 


ln(^ 

ao 


f{N,ao)=MN). 


Then the recursion relations become 


4+i( 7V) = - p(o)(iv)4(iV) = Ro{N)An{N), 

Po 


(165) 


(166) 


(167) 


and can be solved as 


An — 


Ro{N) Ao{N) 


= (e+rl+e-rl) f{N,ao), 


( 168 ) 


having used eq. (nni. Inserting this expression into eq. (116511 we easily obtain the relations 


f{N, a,) = Y. 


[fi„(iV) 


n=0 


n\ 


ln(^ 

ao 


^(iV)= (^j /(iV,«o), 


(169) 


in agreement with eq. (Tmnii . 
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9.2 The standard NLO solution from moment space 


Moving to NLO, one can build a truncated solution in moment space of eq. (Unni) by a 
series expansion around the lowest order solution. 

We start from the truncated version of the vector equation & 


df{N, as) 


das PoOis 

that we re-express in the form 


p(0) < ^ f p(i) _ ^p(0) 

27r V 2 


f{N,as), 
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(y.g 


pm R„b, 


m,a.) 


Rq c^gRi 


f{N,as), 


and with the Ri operator defined as 




(170) 


(171) 


( 172 ) 


We use dacni a truncated vector solution of (nni) - accurate at 0{as) - of the form 

f{N,as) = U{as,N) L{as,ao,N)U-\ao,N) f{N,ao) 


1-1-Q(sf/i(iV) L{as,ao,N) 1-|-Q(oL^i(W) /(W, ckq), 


-1 


(173) 


where we have expanded in powers of as the operators 17(q(s, N) and IJ ^(ckq, N). Inserting 
dnsi) in eq. (dH), we obtain the commutation relations involving the operators t/i, Rq 
and .Ri 


R^,Ui 


IR-R 


1 ) 


(174) 


which appear in the solution in Mellin space dacEi. Then, using the properties of the 
projection operators 


tj\ — e_|_t/ie_|_ -f- C—UiC—, 

and inserting this relation in the commutator di we easily derive the relation 

e+i?ie_ e_i?ie+ 


( 176 ) 


Ui = 


e+i?ie+ -t- e-RiC- 


r+ — r_ 


r_ — 


(176) 
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Finally, expanding the term 


1 + CXqU 1 


in eq. (Hisi) we arrive at the solution US EE] 


An, a,) = 


L -f- oigtj\L — cxoLIJi 


f{N,ao), 


(177) 


where the (a^, Oq, N) dependence has been dropped. Such solution can be put in a more 
readable form as 

^+ 


fiN,ag) = 


a. 


tto 


(Xg 

CXO - CXs \ — 
ao 


e_|_ + (cTs — CTo) e_|_i7ie_|_ + 

e+i?ie_ 
r+ — r_ — 1 


r _— r + 


+ (+<-»•—)[■ /(iV, ao), (178) 


which can be called the standard NLO solution, having been introduced in the literature 
about 20 years ago ng. It is obvious that this solution is a (first) truncated solution of 
the NLO singlet DGLAP equation, with the equation truncated at the same order. 


9.3 Reobtaining the standard NLO solntion using the logarithmic 
expansion 


Having worked out the well-known NLO singlet solution in moment space, our aim is to 
show that the same solution can be reconstructed using a logarithmic ansatz. This fills a 
gap in the previous literature on this types of ansatze for the QCD pdf’s. To facilitate our 
duty, we stress once more that the type of recursion relations obtained in the non-singlet 
and singlet cases are similar. In fact the matrix structure of the equations doesn’t play 
any role in the derivation due to the linearity of the ansatz in the (vector) coefficient 
functions that appear in it. 


Our NLO singlet ansatz has the form 




NLO 


oo 


E 

n=0 



In 





(179) 


with An and Bn now being vectors involving the singlet components. The recursion 
relations are 


An+,{N) = Ro{N)An{N), 

Bn+l{N) = -Bn{N) - ^An+l{N) + Ro{N)Bn{N) - A^pi^){N)An{N), (180) 

47r Trpo 

subjected to the initial condition 

An, ao) = A,{N) + ao5o(7V). (181) 

The solution of dnni) in moment space can be easily found and is given by 

4+i(JV) = [e+(r+)” + e_(r_)”] 4(iV), (182) 
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while we can re-arrange the Bn+i{N) relation into the form 


Bn+i{N) = (^0 - l) Bn{N) + RiR^AoiN). (183) 


The last step to follow in order to construct the truncated solution involves a projection 
of the recursion relations (HHi in the basis of the projectors e±. In this basis we separate 
the equations as 


Bn+i = (e+r+ e_r_ - 1) 


B++ + -h B' 

I -^r) I n 


?+- 


+ 


Rt+ + Rf- + + Rr {e+rl + e.r^) Iq, (184) 


where we have used the notation 


e±Bn{N)e± = B^^. ( 185 ) 

Then, using 

B**, = (r+ - 1) B++ + Rt*rlA„, 

= (r+ - 1) + Rt-r1A„, 

4+t = (r- - 1) B-+ + R^rlAa, 

Bn+i = (’■- - 1) B~~ + R^~r1Aa, (186) 

it is an easy task to verify that the solutions of the recursion relations at NLO are given 
by 

Bi* = [tI - (r+ - 1)”] fl++4, 

B-- = [r" - (r- - 1)"] fir-Ao, 

B*- = [-r^ + (r+ - 1)”] 

B-+ = [-rl + (r_ - 1)"]- 

(187) 


where we have expressed the nth iterate in terms of the initial conditions, and we have 
taken Bq = 0. Summing over all the projections, we arrive at the following expression for 
the NLO truncated solution of the singlet parton distributions 


m a.) = 


n=0 


nl 


Ai* + A-- + a. ( B++ + B-- + B-+ + Bl 


(188) 


which can be easily exponentiated to give 

/(iV, Cts) = C+Aq f -^ -|- e_y4o (—^ + 

\Q( 0 / V'^ 0 / 
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(189) 


Finally, organizing the varions pieces we obtain exactly the solntion in eq. dllHI). We 
have therefore shown that the logarithmic ansatz coincides with the solution of the singlet 
DGLAP equation at NLO known from the previous literature and reported in the previous 
section. It is intuitively obvious that we can build with this approach truncated solutions 
of higher orders improving on the standard solution (HiHi) known from moment space, 
and we can do this with any accuracy. However, before discussing this point in one of the 
following sections, we want to show how the same strategy works at NNLO. 


9.4 Truncated Solution at NNLO 


At this point, to complete our investigation, we need to discuss the generalization of 
the procedure illustrated above to the NNLO case. As usual, we start from a truncated 
version of eq. (nsni), that at NNLO can be written as 


df{N, as) 

das 


1 

Otg 


Ro + agRi + a^R2 


nN,as), 


(190) 


where 


R9 = - 


1 ( P(2) 


+ 


+ 


lOvr 


(191) 


TT \ 2TTj3o 4 
whose solution is expected to be of the form PI 

f(N,a,) = [l + a,i7i(]V) + aJC/2(]V)l L(a„a„,N) \l + a„Ui(N) + alU2(N)Y' 


(192) 


where 


Ro, IJi 

Ro, U2 


= Ui- Pi, 

= —R 2 — RiUi + 2U2- 


(193) 
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Using the projectors in the ± subspaces, one can remove the commutators, obtaining 


= - 
2 2 


U2~ = - 
^ 2 


= 


D++ D++ I D++ 


+ R 2 

1 


r_ — r+ — 1 

R^+Rf- 

r+ — r_ — 1 




r+ — r_ — 2 
1 

r_ — r+ — 2 




- RJ~ + 


-R^^Rt^ - R-, 


r+ — r_ — 1 

R-R^-^ 

r_ — r+ — 1 


and the formal solution from Mellin space can be simplified to 


(194) 


/(TV, as) =[L + asUd - aoLlR 

WsdL - + alL (ul - Ua)! /(iV, ao) • (195) 


At this point we introduce our (1-sf truncated) logarithmic ansatz that is expected to 
reproduce (119511 . Now it includes also an infinite set of new coefficients U„, similar to the 
non-singlet NNLO case 


/A.a.)=5T 

n=0 



-^n ~t“ “t“ CUgCji 


(196) 


Inserting the NNLO logarithmic ansatz into (i™ . we obtain in moment space the 
recursion relations 


Rn+l RoRni 

Bn+1 = (^Ro ~ 1^ Bn + RiRqAq, 


Cn+l — 


Ro — 2) Cn — ( Bn + Bn+l ) + 


-— Rq + Ri 
dvr 


Br, 


R2 + ^Ri 

dvr 




(197) 


whose solution has to coincide with (II9511 . Also in this case, as before, we use the e± 
projectors and notice that the structure of the recursion relations for the coefficients An 
and Bn remain the same as in NLO. Therefore, the solutions of the recursion relations 
for An and Bn are still given by (HHi and HHII). We then have to find only an explicit 
solution of the relations for the new coefficients C'„+i(A^). 

These relations can be solved in terms of Aq, Bq and Cq with the help of (UEZD- Finally, 
taking Bq = 0 and Uq = 0 after a lengthy computation we obtain the explicit solutions 
for the projected components 


a 


++ 


1 _ Rt~Rd _ 

2 (r+ — r_ — l)(r_ — r_|_ — 1) 
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[2(r_ - 1)" - (r+ - 2)” - (r+ - 2)"r+ - rj + r”+‘ + r_ ((r+ - 2)" - r")] A, 
+\Rt*Rr [>■: - 2(r+ - 1)” + (r+ - 2)”] A„ 

K - (>■+ - 2 )”] 4, 


CZ- = -i 




X 


2 (r+ — r_ — l)(r_ — r+ — 1) 

[2(r+ - - (r_ - 2)" - (r_ - 2)”r_ - + rl+^ + r+ ((r_ - 2)^ - r")] Iq 

[r- - 2(r_ - 1)” + (r_ - 2)”] A„ 

+ 54 - - (>■- - 2 )"] 4. 


C.+- = 


2 + + r_(3 — 2r+) — 3r+ + r 


(r” — (r_|_ — 2)”) (1 + r_ — r_|_) 


+ 


-Rf-Rr (2(r- - 1)” + (r_ - 1)V_ - 

-(r+ - 2)*^ + r"(r+ - 1) - (r_ - 1)V+) 

+R++R+- (r" + (r+ - 2)^ + r_(r+ - 2)” 

-2(r+ - 1)” - r_(r+ - I)*" - r+(r+ - 2)” + r+(r+ - I)*")] Aq, 

(?-+ =- 5 - ^5 -^-- fe+(r_ - r+ - 1) ((r_ - 2)^ - r”) 

2 + r2 +3r+ + r^-r_(3 + 2r+) L ^ ^ ^ +) 

(-(r- - 2)" + 2(r_ - 1)” + ((r_ - 2)" - (r_ - l)’")r_ 

-(r_ - 2)V+ + (r_ - 1)V+ - r") 

+Rr+R^+ ((r_ - 2Y - 2(r+ - + r_(r+ - 1)” 

-(r+ - l)"'r+ + r” - r-r'^ + Aq. 


(198) 


Re-inserting these solntions into the NNLO ansatz (IT^ and after exponentiation one can 
show explicitely that the logarithmic solution so obtained coincides with (I195jl . Details 
can be found in appendix A. In the practical implementations of these solutions, there are 
two obvious strategies that can be followed. One consists in the implementation of the 
recursion relations as we have done in various cases above: a sufficiently large number of 
iterates will converge to the truncated solution (If 951) . This is obtained by implementing 
eqs. dMl) and incorporating them into (ITTIHl) . A second method consists in the direct 
computation of {MI) in x-space which becomes 


f{x,as) = 


L + agUi ® L — aoL ® Ui 

-\-oi^U2 ® L — otgCXQUi ® L ® Ui + ch'qL ® ® Ui — U 2 


f{x,ao), 

(199) 


and with L now replaced by its operatorial (®) form 


I = 



( 200 ) 
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An implementation of (ITO would reduce its numerical evaluation to that of a sequence 
of LO solutions built around “artificial” initial conditions given by f/i ® f{x,ao), Ui ® 
Ui ® f{x, Q!o) and so on. 


10 Higher order logarithmic approximation of the NNLO 
singlet solution 

The procedure studied in the previous section can be generalized and applied to obtain 
solutions that retain higher order logarithmic contributions in the NLO/NNLO singlet 
cases. The same procedure is also the one that has been implemented in all the existing 
codes for the singlet: one has to truncate the equation and then try to reach the exact 
solution by a sufficiently high number of iterates. On the other hand, x-space (non brute 
force) implementations are, from this respect, still lagging since they are only based on the 
Rossi-Storrow formulation which we have analized thoroughly and largely extended 
in this work. 

Therefore, the only way at our disposal to reach from x-space the exact solution is by 
using higher order truncates. This is of practical relevance since our algorithm allows to 
perform separate checks between truncated solutions of arbitrary high orders built either 
from Mellin or from x-space We recall, if not obvious, that in the analysis of hadronic 
processes the two criteria of using either truncated or exact solutions are both acceptable. 

To summarize: since in the singlet case is not possible to write down a solution of the 
DGLAP equation in a closed exponential form because of the non-commutativity of the 
operators Ri, the best thing we can do is to arrange the singlet DGLAP equation in the 
truncated form, as in eqs. (nnD and (fimi . The truncated vector solutions (Hini) and (CUD 
are equivalent to those obtained using the vector recursion relations at NLO/NNLO. 

Now, working at NNLO, we will show how the basic NNLO solution can be improved 
and the higher truncates identified. Clearly, it is important to show explicitely that these 
truncates, generated after solving the recursion relations, can be rewritten exactly in the 
form previously known from Mellin space. We are going to show here that this is in fact 
the case, although some of the explicit expressions for the higher order coefficient functions 
Cn, D/s will be given explicitely only in part. The expressions are in fact slightly lengthy 

Therefore, here we will just outline the procedure and illustrate the proof up to the 
second truncate of the NNLO singlet only for the sake of clarity. 

®The current benchmarks available at NNLO are limited to exact solutions and do not involve com¬ 
parisons between truncated solutions 

^They will be included in a file that will be made available in the same distribution of our code, which 
is in preparation. 
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The exact singlet NNLO equation in Mellin space is given by 


dmas) _ (t) PW(iV) + (t)^P«(iV) + 


da. 




f{N,a, 


P 


NNLO 


{N,as 


(3NNLO(^Cis 


-ma,), 


(201) 


where we have introduced the singlet kernels. After a Taylor expansion of jpNNLOf^^^'^ 

up to af it becomes 


df{N,a,) 1 

da, a. 


Rn 


1)2 ^ '' 2 '' ^1 
—— :^Ri + agRi + cr^i?2 — z — agR 2 
[AtiY dvr 




( 202 ) 


which is the truncated equation of order The i?j (i = 0,1, 2) operators are listed below 




R2 = 


TlPo 271^0 
P( 2 ) 


+ 


2n‘^Po Ati^Pq 


.p(i) + 


bl 


+ 


Stt'^Po Stt'^Pq 


p(0) 


The formal solution of this equation can be written as mm 


(203) 


f{N,aY = r„ 


exp 


da' — I R 


' ao 


^ O''. 


^2 fS , _ / A , _ / 2 7^ /3 


(dvr; 


—i?! + a'gR\ + CTg R 2 — TZ.'^'s -^2 


dvr 


= f/(iV, a,)L(a„ ao)f/"niV, «o)/(iV, ^o), 


/(iV,«o) 

(20d) 


where the operator acts on the exponential similarly to a time-ordered product, but 
this time in the space of the couplings. Again, expanding the U{N,as) and U~^{N,ao) 
operators in the formal solution around = 0 and oq = 0 , we have 


/ {N, cTs) — L P agUiL — aoLUi P a^U2L — asaoUiLUi -1- a^L (— U 2 


PaJJ^L + aga^UiL (^P]^ — U2j — 

-all (Uf - U,U 2 - U 2 U 1 P Pa)! f{N,ao). 


(205) 


Inserting the expanded solution into eq. (i™ and equating the various power of a, we 
arrive at the following chain of commutation relations 


Rq, Ui 

Rq, U 2 
Ro, U 3 


Ui - Ri, 

—R 2 — PiPi + 2 P 2 , 

71—71—7“^2 ~ R1U2 — R2U1 + 3P3. 
(dvr)^ (dvr) 


(206) 


d3 
































Removing the commutators by the e± projectors one obtains 


Ut^ = , 

ur = 4", 

ut- = 


Rt 


r+ — r_ — 1 ’ 
r_ — r I — 1 ’ 




dH-+ dH— h I dH— h 


t/o— = ^ 


ut = 


R]^ R^ - 1 “ R 2 


Rj-R^^ 

r_ — r+ — 1 

R^+Rt~ 

r — r_ — I 


ur = 


r+ — r_ — 2 
1 

r_ — r+ — 2 
1 


-Rt~k[~ - Rt~ + 
-R{^Rt^ - + 


Rt"^Rl~ 

r_|_ — r_ — 1 

R~k[+ 

r_ — r+ — 1 


r/++ — 

' “ 3 L (47r) 


^2 


+ Rt~U2^ + + Rt^ut^ 


(47j-'i2 


1 


“ 3 L (47r) 


hi A__ 62 


^r" - 7j^R2~ + + Ri~U2~ + R2^ut~ + R^Ui 


(47r 
hi 


^ r+ — r_ — 3 [ ( 47 r) 

IJ-+ = _ \ __^ 

^ r_ — r+ — 3 ( 47 r) 


&2 


Rt~ - vj^Rt - Rt~ur - Rt^ut - Rt~ur - Rt^ut 


RR 


(4vr) 

h2 

(4vr) 


, 4 -+ - R^+u++ - Rru2^ - R2^ut+ - Rrui 


(207) 


As one can see, the t /3 operator is expressed in terms of the kernels p(o),pd) and P^‘^\ 
One can prove that by imposing a higher order ansatz in Mellin space of the form 


~~ r 

/(TV, a.) = V ^ 4(iV) + asB^N) + a^d(TV) + «'^n(iV) 
^ n! L 


(208) 


n=0 


the solution (IM is generated. The vector recursion relations in this case become 

4+1 = -^p®4. 


Po 


Bn+l — —Bn 


Cn+l — 


X _ T_ 

n ^ -^n n ^n+l n ^ ^rii 

-KfJo TT/To fJo 


2 vr2/To 

A 


p(2)4 


4+1 - 


(4vr)2/3o 




Du-\-i 


_ \ _ — _ B — 

27r2/To " " 


TT/^o' 

/- 

^2 


(dvr) 2/To " (dvr) 2/To tt^o 


- ^p(^)d 
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(209) 


A 

27r/3o 


a 


(47r)/3o Po 


Applying the properties of the e± operators, we project out the ... components 

of these relations, and imposing the initial conditions Bq = Cq = Dq = Q together with 


fiN,ao)=Ao, ( 210 ) 

we obtain the explicit form of /(iV, a^). In order to construct the /(iV, a^) solution, 
the four ± projections of Dn+i must be solved with respect to Ao(iV), since the other 
projections A^, are already known. A direct computation shows that the structure 

of the solution can be organized as follows in terms of the components of Ri 


D++(7V) = r+, r_, N, 4) + W^iRf-Ri^Rf^, r+, r_, N, 4) 

+Ws{Rt~RrRl^, r+, r-, N, 4) + W4(^4^4^r^, r+, r-, N, 4) 
+W,{Rt+Rt^, r+, r_, N, 4) + W,{Rt^Rt+,r+, r_, N, 4) 
+W7(^44-+, r+, r_, N, 4) + WsiRtRl^, r+, r., N, 4) 
+W9(^4, r+, r_, N, 4) + Wio(4++, ^+, r., N, 4), 


D--{N) = Wi(4--3, r+, r_, N, 4) + W^iR^^RtRr ,r+, r-, N, 4) 

+W3(4+^4^4, r+, r_, N, 4) + W,{R~R^+Ri-,r+, r_, N, 4) 
+W 5 (.R]“".R^“, r+, r_, TV, 4) + WeiR 2 ~ Ri~ ,r+, r_, TV, 4) 
+W7(^r+^4, r+, r_, N, 4) + W8(4-+^4, r+, r_, N, 4) 

+W9(.R^ , r_|_, r_, TV, 4) + Wio(4 A+A-:R^:Aq), 


D^-{N) = Zi(^4^4^r-,r+,r_,iV,4) + Z^iRt^Rf^Rt ,r+,r_, N, Ao) 

+Z3(.Rj'“ 4 -^1 Aq) + Z4{R'^ 4 "'"-Ri" 4) 

+Z 5 {R^~R 2 ~,r+, r_, N, 4) + Z(i{R^~ R];~ ,r+, r_, N, 4) 

+Z7(^4^4, r+, r_, iV, 4) + Z8(^4^4, r+, r_, iV, 4) 

+Z9(.Rj'' , r+, r_, iV, 4) + Zio(.R^ , r+, r_, iV, 4); 

D-+{N) = Zi(4-4-+^4, r+, r_, N, 4) + Z^iRr^'Rl^, r+, r., N, 4) 
+Z3(4+i?4^4, r+, r_, N, 4) + Z4(4+^44+, r+, r_, N, 4) 
+Z5(4+i?4, r+, r_, N, 4) + Z6(4”^^4, r+, r_, N, 4) 

+Z7(4"^^’^, r+, r-, N, 4) + Z8(4~^r’^> ^o) 

+Z^{Ri r+, r_, iV, 4) + Zio(4 > ^o); 

( 211 ) 

where we used the notation The expressions of the functions W 
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and Z’s can be extracted by a symbolic manipulation of the coefficients . We have 
included one of the projections for completeness in an appendix for the interested reader. 

In x-space, the structure of the (2nd) truncated (or NNLO solution can be ex¬ 
pressed as a sequence of convolution products of the form 


Os) = L{x) -f asUi{x) 0 L{x) - aoL{x) 0 Ui{x) + alU 2 {x) ® L{x 

—asaoUi{x) ® L{x) ® Ui{x) -f aQL{x) ® (Ui{x) ® Ui{x) — U 2 {x) 

+allJ^{x) ® L{x) + asaltJi{x) ® L{x) ® ^f/i(a;) ® lJi{x) - U 2 {x) 

—alaoU 2 {x) ® L{x) ® Ui{x) — alL{x)® 

Ui{x) ® lJi{x) ® lJi{x) - lJi{x) ® U 2 {x) - U 2 {x) ® Ui{x) + lJ‘i{x 


fix,ao), 

( 212 ) 


and is reproduced by the af logarithmic expansion 


f{x, as) =y^^—r An{x) + asBn{x) + a‘lCn{x) + aliDn{x) 
n\ L 


n=0 


with the initial condition 


(213) 


f(x,ao) = Ao{x). 


(214) 


The study of higher order truncates is performed numerically with the implementation of 
the generalized recursion relations given in the previous sections. 


11 Comparison with existing programs 

In this section we present a numerical test of our solution algorithm. For this aim, we 
compare the results of a computer program that implements our method with the results 
of QCD-Pegasus j2j, a PDF evolution program based on Mellin-space inversion, which 
has been used by the QCD Working Group to set some benchmark results [201 El]. In 
the following, we refer to these results as to the benchmark. 

In the tables we are going to show in this paper, we set the renormalization and 
factorization scales to be equal, and we adopt the fixed flavor number scheme. The final 
evolution scale is /i|. = lO'^GeV^. We limit ourselves to this case because it is enough to 
test the reliability of our method. A more lenghty and detailed analysis, which will take 
into account many other cases with renormalization scale dependence and variable flavor 
number scheme, will be presented elsewhere. 

As in the published benchmarks, we start the evolution at = 2GeV^, where the 
test input distributions, regardless of the perturbative order, are parametrized by the 
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following toy model 


XUy(x) - 
xdy{x) - 

xg{x) ■ 
xd{x) - 
xu{x) ■ 
xs{x) = xs{x) - 

and the rnnning conpling has the valne 


5.107200a;°-®(l - 
3.064320a;°-®(l - xf 
1.700000a;-°-^(l -x)® 

0.1939875x“°-^(l -x)® 

(1 — x)xd(x) 

0.2x{u + d){x) (215) 


asifJ.%0 = 2 GeV^) = 0.35. (216) 

We remind that = q — q, q+ = q + q, L± = d ± u. Onr resnlts are obtained nsing 
the exact solution method for the nonsinglet and the LO singlet, and the K-th truncate 
method for the NLO and NNLO singlet, with k = 10. In each entry in the tables, the 
first number is our result and the second is the difference between our results and the 
benchmark. 

In Table El we compare our results at leading order with the results reported in Table 
2 of (20]. The agreement is excellent for any value of x for the nonsinglet {xUy and xdy ); 
regarding the singlet {xg column), the agreement is excellent except at very high x: we 
have a sizeable difference at x = 0.9. In Table IH we analyze the next-to-leading order 
evolution; the results for the proposed benchmarks are reported in Table 3 of (SHI- The 
agreement is very good for any value of x for the nonsinglet; for the singlet the agreement 
is good, except at very high x (x = 0.9). 

Moving to the NNLO case (Table El the benchmarks are shown in Table 14 of ini). 
some comments are in order. We don’t solve the non-singlet equation as in PEGASUS, 
since (inn)) admits an exact solution jnHD, which in PEGASUS is obtained only by itera¬ 
tion of truncated solutions. Our implementation is based on the exact solution presented 
in this work. The discrepancy between our results and PEGASUS are of the order of 
few percent, and they become large in the gluon case at x = 0.9, as for the lower or¬ 
ders. Another comments should be made for the sea asimmetry of the s quark, which is 
nonvanishing at NNLO. In this case we have a sizeable relative discrepancy for any value 
of X (reported in the column xs„), but it is evident that this asymmetry is quite small, 
especially if compared with the column xs+, whose entries are several orders of magnitude 
larger. This means that xs and xs should be comparable and their difference therefore 
pretty small. We have to notice that xsy cannot be computed directly by a single evo¬ 
lution equation. Indeed, it is computed by a difference between very close numbers, a 
procedure that amplifies the relative error. 
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LO 

, Uf = 

4, = 10^ GeV2 






X 

XUi 

) 

xdy 

xL. 


2xL 

+ 

XS+ 

XC+ 

xg 

10“^ 

5.7722 • 

10-6 

3.4343 • 

10-6 

7.6527 • 

10-^ 

9.9465 • 

10+6 

4.8642 • 

10+6 

4.7914 • 

10+6 

1.3162 • 

10+6 


0.0000 • 

10-6 

0.0000 • 

10-6 

0.0000 • 

10-7 

0.0000 • 

10+6 

0.0000 • 

10+6 

0.0000 • 

10+6 

0.0000 • 

10+6 

10-6 

3.3373 ■ 

10-“^ 

1.9800 • 

10-"^ 

5.0137 • 

10-6 

5.0259 ■ 

10+6 

2.4263 • 

10+6 

2.3685 ■ 

10+6 

6.0008 • 

10+2 


0.0000 • 

10--^ 

0.0000 • 

10-^ 

0.0000 • 

10-6 

0.0000 • 

10+6 

0.0000 • 

10+6 

0.0000 • 

10+6 

0.0000 • 

10+2 

10-6 

1.8724 • 

10-6 

1.1065 • 

10-6 

3.1696 • 

10-6 

2.4378 • 

10+6 

1.1501 • 

10+6 

1.1042 • 

10+6 

2.5419 • 

10+2 


0.0000 • 

10-6 

0.0000 • 

10-6 

0.0000 • 

10-6 

0.0000 • 

10+6 

0.0000 • 

10+6 

0.0000 • 

10+6 

0.0000 • 

10+2 

10 -^ 

1.0057 ■ 

10-2 

5.9076 • 

10-6 

1.9071 • 

10 -^ 

1.1323 ■ 

10+6 

5.1164 • 

10 +° 

4.7530 ■ 

10 +° 

9.7371 • 

10+6 


0.0000 ■ 

10-2 

0.0000 • 

10-6 

0.0000 • 

10 -^ 

0.0000 ■ 

10+6 

0.0000 • 

10 +° 

0.0000 ■ 

10 +° 

0.0000 • 

10+6 

10-6 

5.0392 • 

10-2 

2.9296 • 

10-2 

1.0618 • 

10-6 

5.0324 • 

10+6 

2.0918 • 

10 +° 

1.8089 • 

10 +° 

3.2078 • 

10+6 


0.0000 • 

10-2 

0.0000 • 

10-2 

0.0000 • 

10-6 

0.0000 • 

10+6 

0.0000 • 

10 +° 

0.0000 • 

10 +° 

0.0000 • 

10+6 

10-2 

2.1955 ■ 

10-6 

1.2433 • 

10-6 

4.9731 • 

10-6 

2.0433 ■ 

10+6 

7.2814 • 

10-6 

5.3247 ■ 

10-6 

8.0546 • 

10 +° 


0.0000 ■ 

10-6 

0.0000 • 

10-6 

0.0000 • 

10-6 

0.0000 ■ 

10+6 

0.0000 • 

10-6 

0.0000 ■ 

10-6 

0.0000 • 

10 +° 

0.1 

5.7267 • 

10-6 

2.8413 • 

10-6 

1.0470 • 

10-2 

4.0832 • 

10-6 

1.1698 • 

10-6 

5.8864 • 

10-2 

8.8766 • 

10-6 


0.0000 • 

10-6 

0.0000 • 

10-6 

0.0000 • 

10-2 

0.0000 • 

10-6 

0.0000 • 

10-6 

0.0000 • 

10-2 

0.0000 • 

10-6 

0.3 

3.7925 • 

10-6 

1.4186 • 

10-6 

3.3029 • 

10-6 

4.0165 • 

10-2 

1.0516 • 

10-2 

4.1380 • 

10-6 

8.2676 • 

10-2 


0.0000 ■ 

10-6 

0.0000 • 

10-6 

0.0000 • 

10-6 

0.0000 ■ 

10-2 

0.0000 • 

10-2 

+ 0.0001 

• 10-6 

0.0000 • 

10-2 

0.5 

1.3476 ■ 

10-6 

3.5364 • 

10-2 

4.2815 • 

10 -^ 

2.8624 ■ 

10-6 

7.3137 • 

10 -"^ 

2.6481 ■ 

10-'6 

7.9242 • 

10-6 


0.0000 • 

10-6 

0.0000 • 

10-2 

0.0000 • 

10 -^ 

0.0000 • 

10-6 

- 0.0001 

• 10 -^ 

0.0000 • 

10-'6 

+ 0.0002 

• 10-6 

0.7 

2.3123 • 

10-2 

3.5943 • 

10-6 

1.5868 • 

10-6 

6.8970 • 

10-6 

1.7730 • 

10-6 

6.5593 • 

10 -° 

3.7301 • 

10-6 


0.0000 • 

10-2 

0.0000 • 

10-6 

0.0000 • 

10-6 

+ 0.0009 

• 10-6 

+ 0.0005 

• 10-6 

+ 0.0044 

• 10 -° 

- 0.0010 

• 10-6 

0.9 

4.3443 ■ 

10-‘6 

2.2287 • 

10-6 

1.1042 • 

10-6 

3.3030 ■ 

10-6 

8.5607 • 

10 -° 

3.2577 ■ 

10 -° 

1.3887 • 

10 -° 


0.0000 • 

10-^6 

0.0000 • 

10-6 

0.0000 • 

10-6 

- 0.3263 

■ 10-6 

- 0.1631 

• 10-6 

- 1.6316 

• 10 -° 

+ 0.2969 

• 10 -° 


Table 3: Comparison between our and the benchmark results at LO. In each entry, the first number is our result and the second one is 
the difference between our result and the benchmark. The benchmark values are reported in Table 2 of 1201 































NLO, Uf = 4, = 10^ GeV^ 

X 

XUy 

Xdy 


2xL+ 

XS+ 

XC + 

xg 

10-7 

1.0620 • 10-“^ 
+0.0004 • 10-^ 

6.2353 • 10-6 
+0.0025 • 10-6 

4.2455 • 10-6 
+0.0015 • 10-6 

1.3710 • 10+2 
+0.0112 • 10+2 

6.7469 • 10+6 
+0.0556 • 10+6 

6.6750 • 10+6 
+0.0555 • 10+6 

1.1517 • 10+6 
+0.0034 • 10+6 

10-6 

5.4196 • 10-^ 
+0.0019 • 10-^ 

3.1730 • 10-^6 
+0.0011 • 10-^ 

1.9247 • 10-6 
+0.0006 ■ 10-6 

6.8896 • 10+6 
+0.0500 • 10+6 

3.3592 ■ 10+6 
+0.0250 • 10+6 

3.3021 • 10+6 
+0.0250 • 10+6 

5.4048 ■ 10+2 
+0.0137 • 10+2 

10-6 

2.6878 • 10-6 
+0.0008 • 10-6 

1.5682 • 10-6 
+0.0005 • 10-6 

8.3598 • 10-6 
+0.0023 • 10-6 

3.2936 • 10+6 
+0.0208 • 10+6 

1.5788 • 10+6 
+0.0103 • 10+6 

1.5335 • 10+6 
+0.0104 • 10+6 

2.3578 • 10+2 
+0.0050 • 10+2 

10-^ 

1.2844 • 10-2 
+0.0003 • 10-2 

7.4576 • 10-6 
+0.0018 • 10-6 

3.4919 • 10-"^ 
+0.0008 ■ 10-^ 

1.4824 • 10+6 
+0.0078 • 10+6 

6.8744 • lO+o 
+0.0389 • 10+9 

6.5156 • 10+9 
+0.0387 • 10+9 

9.3026 • 10+6 
+0.0154 • 10+6 

10-6 

5.7937 • 10-2 
+0.0011 • 10-2 

3.3343 • 10-2 
+0.0006 • 10-2 

1.4164 • 10-6 
+0.0002 • 10-6 

6.1899 • lO+o 
+0.0251 • 10+9 

2.6783 ■ lO+o 
+0.0124 • 10+9 

2.4001 • 10+9 
+0.0123 • 10+9 

3.1540 ■ 10+6 
+0.0038 • 10+6 

10-2 

2.3029 • 10-6 
+ 0.0003 • 10-6 

1.2930 • 10-6 
+ 0.0002 • 10-6 

5.3258 • 10-6 
+ 0.0007 • 10-6 

2.2587 • lO+o 
+ 0.0060 • 10+9 

8.4518 • 10-6 
+ 0.0298 • 10-6 

6.5540 • 10-6 
+ 0.0294 • 10-6 

8.1120 • 10+9 
+ 0.0054 • 10+9 

0.1 

5.5456 • 10-6 
+ 0.0004 • 10-6 

2.7338 • 10-6 
+ 0.0002 • 10-6 

1.0012 • 10-2 
+ 0.0001 ■ 10-2 

3.9392 • 10-6 
+ 0.0056 • 10-6 

1.1517 ■ 10-6 
+ 0.0028 • 10-6 

6.0619 • 10-2 
+ 0.0268 • 10-2 

8.9872 ■ 10-6 
+ 0.0005 • 10-6 

0.3 

3.5395 • 10-6 
+ 0.0002 • 10-6 

1.3158 • 10-6 
0.0000 • 10-6 

3.0363 • 10-6 
+ 0.0001 • 10-6 

3.5884 • 10-2 
+ 0.0036 • 10-2 

9.2210 • 10-6 
+ 0.0180 • 10-6 

3.4066 • 10-6 
+ 0.0176 • 10-6 

8.3415 • 10-2 
- 0.0036 • 10-2 

0.5 

1.2271 • 10-6 

0.0000 • 10-6 

3.1968 • 10-2 
+ 0.0001 • 10-2 

3.8266 • 10 -^ 
+ 0.0001 ■ 10 -^ 

2.4149 • 10-6 
+ 0.0023 • 10-6 

5.8539 ■ 10 -"^ 
+ 0.0115 • 10 -^ 

1.7068 • 10 -^ 
+ 0.0113 • 10-6 

8.0412 ■ 10-6 
- 0.0061 • 10-6 

0.7 

2.0429 • 10-2 
0.0000 • 10-2 

3.1474 • 10-6 
+ 0.0001 • 10-6 

1.3701 • 10-6 
0.0000 • 10-6 

5.3703 • 10-6 
+ 0.0081 • 10-6 

1.2432 • 10-6 
+ 0.0039 • 10-6 

2.8201 • 10-9 
+ 0.0394 • 10-9 

3.8654 • 10-6 
- 0.0067 • 10-6 

0.9 

3.6097 • 10-6 
+ 0.0001 • 10-^6 

1.8317 • 10-6 
0.0000 • 10-6 

8.9176 • 10-9 
- 0.0054 ■ 10-9 

1.6196 • 10-6 
- 0.4724 • 10-6 

1.6717 • 10-9 
- 2.3673 • 10-9 

- 2.6084 • 10-9 
- 2.3681 • 10-9 

1.8308 • 10-9 
+ 0.6181 • 10-9 


Table 4: Same as in Table 01 in the NLO case. The benchmark valnes are reported in Table 3 of I2DI 
























NNLO, Uf = 4, = 10^ GeV^ 

X 

XUy 

Xdy 


2xL+ 

XSy 

XS + 

XC+ 

xg 

10-7 

1.4069 • 10-^ 
-0.1218 • 10-“^ 

9.0435 • 10-3 
-0.1201 • 10-7 

5.5759 • 10-3 
-0.1259 • 10-3 

1.4184 • 10+2 
+0.0994 • 10+2 

1.9569 • 10-3 
-1.1868 • 10-3 

6.9844 • 10+7 
+0.4967 • 10+7 

6.9127 • 10+7 
+0.4966 • 10+7 

1.0519 • 10+3 
+0.0543 • 10+3 

10-® 

6.5756 • 10-"^ 
-0.3420 • 10-“^ 

4.0826 ■ 10-7 
-0.3458 ■ 10-7 

2.4722 • 10-3 

-0.0688 • 10-3 

7.1794 ■ 10+7 
+0.3295 • 10+7 

5.8862 • 10-3 
-3.5417 • 10-3 

3.5043 • 10+7 
+0.1646 • 10+7 

3.4474 ■ 10+7 
+0.1646 ■ 10+7 

5.1093 • 10+2 
+0.1969 • 10+2 

10-® 

3.0260 • 10-3 
-0.0721 • 10-3 

1.8199 • 10-3 
-0.0775 • 10-3 

1.0393 • 10-7 
-0.0326 • 10-7 

3.4373 • 10+7 
+0.0902 • 10+7 

1.4230 • 10-7 
-0.8560 • 10-7 

1.6509 • 10+7 
+0.0450 • 10+7 

1.6057 • 10+7 
+0.0450 • 10+7 

2.2899 • 10+2 
+0.0602 • 10+2 

10-^ 

1.3656 • 10-2 
-0.0066 • 10-2 

8.0052 • 10-3 
-0.0967 ■ 10-3 

4.1299 • 10-7 
-0.1259 • 10-7 

1.5403 • 10+7 
+0.0199 • 10+7 

2.2837 • 10-7 
-1.3807 • 10-7 

7.1661 • 10+9 
+0.0991 • 10+9 

6.8085 • 10+9 
+0.0988 ■ 10+9 

9.2125 • 10+7 
+0.1457 • 10+7 

10-3 

5.9360 • 10-2 
+0.0200 • 10-2 

3.4135 ■ 10-2 
+0.0085 • 10-2 

1.5650 • 10-3 
-0.0358 • 10-3 

6.3657 ■ 10+9 
+0.0427 • 10+9 

8.9572 • 10-3 
- 0.5522 • 10-7 

2.7684 • 10+9 
+ 0.0210 • 10+9 

2.4913 ■ 10+9 
+ 0.0209 • 10+9 

3.1592 • 10+7 
+ 0.0243 • 10+7 

10-2 

2.3139 • 10-7 
+ 0.0061 • 10-7 

1.2958 • 10-7 
+ 0.0039 • 10-7 

5.5064 • 10-3 
- 0.0624 • 10-3 

2.2868 • 10+9 
+ 0.0116 • 10+9 

- 3.5702 • 10-7 
+ 2.1611 • 10-7 

8.6094 • 10-7 
+ 0.0592 • 10-7 

6.7224 • 10-7 
+ 0.0601 • 10-7 

8.1503 • 10+9 
+ 0.0122 • 10+9 

0.1 

5.5125 • 10-7 
- 0.0052 • 10-7 

2.7142 ■ 10-7 

- 0.0023 ■ 10-7 

9.9834 • 10-3 
- 0.0040 • 10-2 

3.9119 ■ 10-7 
+ 0.0100 • 10-7 

- 1.9045 • 10-7 
+ 1.1582 • 10-7 

1.1453 • 10-7 
+ 0.0067 • 10-7 

6.0520 ■ 10-2 
+ 0.0747 ■ 10-2 

8.9909 • 10-7 
- 0.0654 • 10-7 

0.3 

3.5017 • 10-7 
- 0.0054 • 10-7 

1.3005 • 10-7 
- 0.0020 • 10-7 

3.0025 • 10-3 
- 0.0073 • 10-3 

3.4975 • 10-2 
- 0.0383 • 10-2 

- 1.9830 • 10-3 
+ 1.2061 • 10-3 

8.8758 • 10-3 
- 0.1722 • 10-3 

3.1421 • 10-3 
- 0.1640 • 10-3 

8.3041 • 10-2 
- 0.1145 • 10-2 

0.5 

1.2099 • 10-7 
- 0.0018 • 10-7 

3.1485 ■ 10-2 
- 0.0043 ■ 10-2 

3.7667 • 10-7 
- 0.0075 • 10-7 

2.1876 ■ 10-3 
- 0.1991 • 10-3 

- 1.6924 • 10-9 
+ 1.0291 • 10-9 

4.8155 • 10-7 
- 0.9810 • 10-7 

7.4120 ■ 10-3 
- 0.9758 ■ 10-7 

7.9784 • 10-3 
- 0.1342 • 10-3 

0.7 

2.0052 • 10-2 
- 0.0025 • 10-2 

3.0849 • 10-3 
- 0.0037 • 10-3 

1.3411 • 10-3 
- 0.0023 • 10-3 

1.5984 • 10-3 
- 3.8260 • 10-3 

- 6.2854 • 10-3 
+ 0.3821 • 10-'7 

- 6.1500 • 10-9 
- 1.9084 • 10-3 

- 1.5545 • 10-3 
- 1.9075 • 10-3 

3.8226 • 10-7 
- 0.0722 • 10-7 

0.9 

3.5078 • 10 -^ 
- 0.0033 • 10-7 

1.7767 • 10-3 
- 0.0016 ■ 10-3 

8.6326 • 10-9 
- 0.0184 • 10-9 

- 6.4293 • 10 -'^ 
- 6.6986 • 10 - 7 ' 

- 9.1828 • 10-77 
+ 0.5579 • 10-79 

- 3.2776 • 10 -^ 
- 3.3487 • 10 -'^ 

- 3.3190 • 10 -^ 
- 3.3487 ■ 10 - 7 ' 

1.9280 • 10-9 
+ 0.7144 • 10-9 


Table 5: Same as in Table El in the NNLO case. The benchmark valnes are reported in Table 14 of HQ 

























12 Conclusions 


We have shown that logarithmic expansions identified in x-space and implemented in this 
space carry the same information as the solution of evolution equations in Mellin space. 
This has been obtained by the introduction of new and generalized expansions that we 
expect to be very useful in order to establish benchmarks for the evolution of the pdf’s 
at the LHC. Our analysis has been presented up to NNLO. We have also shown how 
exact expansions can be derived. We have presented analytical proofs of the equivalence 
and clarified the role of previous similar analysis which were quite limited in their reach. 
The numerical implementation of our results will be illustrated in a forthcoming paper 
where various comparisons between our approach and other approaches will be analized 
thoroughly. There are in fact several issues which are still unclear in this area and concern 
the role of the NNLO effects in the evolution and in the hard scatterings, the role of the 
theoretical errors in the determination of the pdf’s, whether they dominate over the NNLO 
effects or not, and the impact of the choices of various truncations in the determination 
of the numerical solution of the pdf’s, along the lines of our work. Similar analysis can 
be performed by other methods, but we think that it is important, in the search for 
precise determination of cross sections at the LHC, to state clearly which algorithm is 
implemented and what accuracy is retained, with a particular attention to the issues 
connected to the resummation of the perturbative expansion 121. Our work, here, has 
been limited to a (fixed order) NNLO analysis. We hope that our analysis has shown 
conclusively that x-space approaches have a very solid base and provide a simple view on 
the structure of the solutions of the DGLAP equations, valid to all orders. We have also 
shown a numerical comparison of our results against those obtained using PEGASUS, for 
a specific setting. The overall agreement, as we have seen, is very good down to very 
small x-values. We will return on this and other related points in future work. 
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13 Appendix A. Derivation of the recnrsion relations 
at NNLO 

As an illustration we have included here a derivation of the recursion relations for the 
first truncated ansatz of O(a^) that appears at NNLO. 

Inserting the NNLO truncated ansatz for the solution into the DGLAP equation we 
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get at the left-hand-side of the defining equation 

OO 

E 


n=l 






n\ as n\ 

2Cn{x) „_]_/3(as) 


a.. 
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-nU 


n\ 


a. 


E \ + 2«,,3(a.)^L" 


n=0 


n\ 


n\ 


( 217 ) 


Note that the first sum starts at n = 1, because the n = 0 term in jnn) does not have 
a dependence. Sending n —> n -|- 1 in the first sum, using the three-loop expansion 
of the beta function 0 and neglecting all the terms of order or higher, the previous 
formula becomes 
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(218) 


At this point we use the NNLO expansion of the kernels. We get at the right-hand-side 
of the defining equation 

OO T C 2 


n=0 


Cto 


[P(2) ^ (^) + ^ [p(0) ^ 


Stt^ 

CK? 


27r 

q; 


dvr^ 


[P'" 0 B„] (x) + E [pO) 0 c„] (x) } . 


(219) 


Equating Jam) and jnn) term by term and grouping the terms proportional respectively 
to as, a^ and a^ we get the three desired recursion relations ii. Setting Q = Qo in jsi 
we get 

f(x, Qo) = Ao(x) + as(Qo)Bo(x) + (as(Q^)f Co(x). (220) 

We have seen that the initial conditions should be chosen as 

Bo(x) = Co(x) = 0, f(x, Qo) = ^o(a;) ( 221 ) 

in order to reproduce the moments of the truncated solution of the DGLAP equation. 


14 Appendix B. NNLO singlet truncated solution. 


Putting all the projections of the coefficients into the NNLO singlet ansatz we get 


OO Y 

f{N , Q(s) = ^ — A++ + A + “s (^B^^ + B^ + + 
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and exponentiating 
f{N, as) = e+Ao 
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This expression is equivalent to that obtained from eq. (nnni). Projecting over all the 
± components we can write 
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Using Ao = f{N, ao) = q{N, ao) we obtain 


a^ 

ao 


r+ 


(226) 


(227) 


g(iV,a,)++ 
q{N,as)~~ 
g(7V,a,)+- 
q{N, as)~^ 


n\ L 


+ alC^ 


n=0 
oo 7- ^ 

E 

n=0 
oo 

E 

n=0 

°° T n ^ 

E TJ ^ 


n! 


a,5+- + a2C'+- 


n=0 


For example we can check the first of the relations above, gives 
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Factorizing (as/ao)'’+ and expanding the power of as the previous expression becomes 
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>++ 


ag 


r+ 


ba + {as - ao)e+Rie+ + -a^ 


e+Rie+Rie-t- + e+i?2e+ 


e+-Rie_i?ie+ 
(r_ — r+ — 1) 


7 O'/-) 


-agO^o 


e+Rie-Rie+ 
(r+ — r_ — 1) 

c-\-Ri&-\-RiC-\- 


+ e+Rie+Rie+ — e+R2e+ 
e+Rie-Rie+ 


a. 


(r+ — r_ — l)(r_ — r+ — 1) \ao 
which agrees with the left hand side of eq. (1^ . 


V— —Vj^ 


f{N,ao). 
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15 Appendix B. Calculation of 

An explicit calculation of the vector coefficient -D++ of the k = A (4th truncated) solution 
of the NNLO singlet equation has been done in this section. Since the expressions of the 
coefficients D~^ and D~~ have a structure similar to we omit them and 

only the explicit form of this one 


give 


D++ = (-3 + r+)" (4 62 - 32 vr^ R++ ^ 


+ 16 5i TT - 5 62 -RC + 40TT^ Ri+^rJ 
-20 61 TT R++ rJ + b 2 Ri^ rJ R^+ ^rJ + Abiir R++ + 96 vr^ i?++ ^ “ 

120 71^ Rf+^rJ -f 247r^ Rf+^r_^ ( 


, - 967r^Aj^+^f l + ^+ 
-3 + r+J ^ \-3 + r^ 

-1 -f r+' 


:++V_ r+ + 


y—j-r/+/ y—j-r/ + / 

120 71^ R^+^r_^ - 24TT^ (^ "3!^^^ ^ + Wb 2 R^^ r_ r+- SOir^ Rf+^ 

40 61 TT R 2 ~^ r_ r+ — 4 62 Ri^ r+ -|- 32 tt^ ^-l- — 16 61 tt r_^ r+-|- 

240 71^ j?++V- r+ - 9671^ R++\J r+- 2407r2 i2++3 r_ r+ 

- 1-96 TT^ r_^ ^ ^ ^ r+ — 5 62 Ri~^ -|- 40 tt^ r+^ — 20 6i tt -|- 662 -RC 

-48 7r^.RC^^-^^+^ + 246i7rR^+r_V+^ - 12071^^5^+^ r+^ 

-H447r^ Rt+^r_^ /^ -2 + ^-r \ 2 ^^20 tt^ R++3 f r+^ - 144 tt^ Rt+^r_^ /^ "l + ^+ \ 2 

V-3-hr+y v-3 + r+y V-3 + r-+y 

—4 62 RC r- r+^ -|- 32 tt^ Ry’*' ^ r_ r+^ — 16 61 tt Rt^ r_ r+^ — 96 tt^ Ry"*" ^ r_ (—^ 3_|_ 

V“3-hr+y 

967r2Ry+^r_ ^ " 3 ^^'*' ^ r+^ -h 62 -Ry^ r+'^ - Stt^ Ry+3 r+^ -h 45i tt Ry+; 

24 fiP » (^ j %/- 24 » (^) ^ 


?y+r+^+ 


^)”,.,‘-4.,i2-(^) 

.40 Rt*^rJ (-^)" + 206. . J (-jT^)” - 62 R+ (-^ 


-f32 7r^Ry+^( /y^ ) - 16 5i TT Ry+f —^ 


.“3 -h r+^ 


+ 


56 




































8 vr^ “ r_ 


(=3T7;) " - 4 0. - 4"* (-^) " - 10 fc r- (^) ’ 


r+ V 


+80 7r2i2++3r_r+- 40 6i vr i?++r_ r ++ Ab 2 Rf^ rJ r+j 


r+ Y 


-32 TT^ i?: 


2 dH— h 3 


‘’■+(=3^7;) +'^'’='^**’'+'(=3T7;)- 

40 .1 j ^ ^44 ” 

+48-71^ r-^ (———^ — 246i tt r_^ r+^ T^ + 4b2 Rt^ r^ (———^ 

V-3 + r+y \-3 + r+J \-3 + r+J 

-32 « r_ 44^ (^) % 16 io . 4. 44=> (^) " - k Rr ^4* ’■ 


+87r^ Rf+'^r+‘^ 


_Y 


-3 + 447 ’■+ 

+16 TT^ (—2 + r_ + r_^ — r_|_ — 2 r_ r+ + r_|_^) x 
-2 + r_\" f r+ \ 


Y 

-3 + r+y 


-2 + r_ + 3 


-87r^ -2 + 3 


-3 + r+y 

2 + r+ Y 


-3 + r+y 

r+ 


— r_ 




r+ 


-3 + r+ 



i?l i?2 


-3 + r+y V-3 + r+_ 

(4 + r_"^ — 4r_^ r+ — 5r+^ + r+"^ + r_^ (—5 + 6r+^) + r_ (lOr+ — 4r+^)) -Ry’*' 
+32 TT^ R+“ R]“+ - 16 TT^ r_ R^" R]“+ - 32 vr^ r.^ R^" R]“+ + 16 tt^ rJ R^~ R]“+ 


-96 TT^ 


) Ry"R^+- 48 7rV_ 

-3 + r+; 2 1 

+16 TT^ r+ R^” R^"*" + 64 vr^ r_ r+ R^” R^"*" — 48 tt^ r_^ r+ R^” Rj^ 


^)"+-+- + 484^4j 


+±^v+-+ 

-3 + r+; 2 1 


+48 TT^ 


— 1 + r_ \ 

_g_^ Y ) - 967rV_ 

—32 TT^ r+^ R^” R]"’*' + 48 vr^ r_ r+^ R^” R^~ 


— 1 + r_ 


-3 + r+y 


Y 


r+ R+ R^ 


+48 77^ 


+64 


44"+-++-16.r24+3ii4-R- 

-^)"+-+-+ 64 .^ 4 R ’■- ^’’ 

-3 + r+y 


-^- I Rt" Rr+ - 16 TT^ rJ 

-3 + r+; 2 1 


r+ Y 

-3 + r+y 


R 2 Ri 


— 16772 r_^ 




R 2 


+32772 r_ r+ 


Y 

-3 + r+y 


Ry R;^ ^ + 48772 r_2 r+ 


r+ Y 
-3 + r+y 


R 2 R]^ 


— 16 772 r+2 




R 2 Ri 


+ 16 772 r+^ 




4 £>++ £>++ 


—96 772 


R++R+++ 120772 r_2 
-3 + r+; 2 1 


R++R++- 24772 r_^ 
-3 + r+; 2 1 


-l + r+\ 


3T^J 


+80 772 r_ r+ R++ R++ - 32 772 rJ r+ Ry+ Ry+ - 240 772 r_ ^ ^2"" 


+96772 ^^3 


r+ Ry+ Ry+ - 40 772 r+2 R++ R++ + 48 772 rJ r+2 R++ R++ 

—3 + r+/ 
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+ 120 TT 

+ 96 TT^ 
-24 TT^ 
-80 vr^ 
+ 160 TT 

-80 TT^ 

— 647 r^ 
+48 TT^ 
-48 TT^ 
-96 TT^ 
+48 TT^ 
32 TT^ r_ 
+ 96 TT^ 
-96 TT^ 
-48 vr^ 
+ 32 TT^ 

+ 16 TT^ 

-32 TT^ 
-32 TT^ 
— 87 r^ r 

-96 TT^ 

+ 24 TT^ 
-32 TT^ 
-48 TT^ 


V-3 + r+y - 

+ 87r2 

IJ__\ n++ 0 ++ 

I -^^2 1 


’■-(5^)“’■4+"+" 

+ 64 ' 

(=3T+) ”+"+" +16 (r5^) ” +" 1^4 

( _3 + ^^ ) Rt^Rf^-M7r^rJr+(^-^^^ Rt^ Rf^ 

4++i?|+ + 96 7r2r_2r+2(^-^^) Rt^ Rf^ 

) n / \ n 

++ ii++ + 16 r+* (z 3 ^) ++ «+ - 32 ir'l +- Jir- -Rf 


’■--'4^33T+j 


. + - 16 TT^ r _2 i ?4 




) Rf-R^--48TT^rJ 
.-3 + r+; ^ 1 1 


Rf kr ++ + 96(r|^) 

4 +^^ R+-R--R- 


;) 


-3 + r+^ 


rJ ( ~ 4 ^~ ) ^+“ i?^+ - 48 TT^ r+ i?^++ 

\—3 + r+y 

-r+^+~ + 487r^ ?•+^+-+ 

’■-( rl ^) ’■++^+^+■1 +48 r ++-+- + 

r_ r++"+“++-16iir+i+■+■++ 

(51$+) ’'4+-+-+++48 ++-+-+- 


(+++)" «4 4- + 48.1 r_ (r$5+)" +- +- 4 

’■4(r$5+) +-+-+*-48 .V+( 35 ^) 


■ + 


Rf~ 


—3 + r+y 

’■-"+ (+$+)” -”4 Rr^ Rr+ +16 .1 r+i (r$5+)" r 4 Rr^ Rf" 

i?$“ + i?$+ + 32 TT^ r_ R+- R^+ R++ + 8 tt^ rJ Rf- + i?$+ 

_ 3 i?+-i?-+i?++ + 967 r 2 (^|^) i?+-i?-+i?++ + 48 7 r 2 r_(^^^) Rf'R^+Rt+ 

1 ^- 1 ^^\ j^+-^-+^++ _ 9 g^ 2 ^_ $ -1 +^+ $ ^+-^ +^++ 

-3 +r+y 1 1 1 V -3 + r+y ^ ^ ^ 

— 1 + r+ \ 

o , „ -^1 Ri Ri 


5^) +24 .ir_4(y^) 11 . 11 . 

■2'■_ V+i?$“ i?y+i?y+ 


'.=4r+i 

r+i?+"i?y+i?++- 167rV_r+i?+“i?j; ,__ 

+ r+i?+-i?r+i?++ 


\ - ’ • -r / 

?$■ i?y+i?y+ + 247r2 r_ 
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-48 7r^r- j r+R^ + i^++ - 72 vr^ r_^ ( -3 + ^ 

+8 TT^ r+2 i?+- R-+ R++ - 24 vr^ r_ r+2 i?+- i?-+ i?++ 

+ r+2 i?+-/?-+i?++ 

+8 TT^ r+3 i?+- i?++ - 24 TT^ ^^3 ^+- ^-+ ^++ ^ 32 ^2 

+16 (z^^) " - 32 vr^ r _2 ( 3 ^^) " Rt' Ri^ R^ 

+64 TT^ r_ r+ ^ ’’+ (^ 3^+ r ) 

-32^'"+'(z^^) ^r^r+^^^-48vr2r_r+2(^-^^) i?r+i?++ 

+16 TT^ r+=^ ^ 3 ^+r ) ~ ^ ^ 


r+2 R+- R-+ R++ 


T^ r+3 R+- R-+ i?++ - 24 TT^ f 

' ' ' V-3 + r+, 




:^+ - 32 TT^ i?++ R+- R-+ + 16 TT^ r_ i?^+ + 

3/?++ii^i?r+- 967r2(^l^) Rt^Rf-R^ 


+32 TT^ r_2 i?J^+ Ri~ R^+ - 16 tt^ rJ i?++ i?+“ i?j“+ - 96 tt^ (^ ^ ^ 

\—3 + r’+y 

+‘‘** ’■- (=FT+) ” R* fl+ -Rr+ + 96(rf^) ”«+ fl+ -Rf^- 
°'^’""’’-( -3 + rl ) +*+■+*-24 +++“+ + 

+24TT^ r_3 ( " 3 ^^^ ) Rf^ Rt~ Ri^ - 16vr^ r+ i?^+ + - 64vr^ r_ r+ R++ R^ 

+48 TT^ rJ r+ R++ R+- R^+ - 48 vr^ r+ Rf+ Rf- R^++ 


>++ pH— f>- 

1 Ri Ri 


967r2 - 

V-3 + r+ 


It+) " ’■* R*R~ R *+(^ 




-72 TT^ r_2 ^ S + r"^ ) 


—48 TT^ r_ r+^ R^ 


i7++j?+~ j?f+- 247r^ (" ^ + ^+ ^ r+2 i?++i?+-i?: 

\—3 + r+y 

r+2i?++7?+-i?-+ +16 7r2r+3 72++i?+-i?r 

3 + r+ / 


^) 7?++7^^7^^ 


+72 TT^ r_ ^ 3 !^^^ ) ^+^ ~ ^i""^ + 16 TT^ r+3 72]^+ 72f+ 

+ 72++ 72+- 72r+ 

+32^'"- (z^^) 72++ 72^72^+ 

72++ 72^72r+ 

+16 7r^r-r+^ ^ ^|+^f+247r^ r+^ ^|+ 

-8^'^+'(z^^) ^^^^r^r+-24vr2r_r+2(^-^^) 72++ 72^ 72^ 


59 



[48 vr^ (4 + r_^ — 4r_^ r+ — 5r+^ + r+^ + r_^ (—5 + 6r+^) + r_ (lOr+ — 4r+^))] ^ . 

(231) 
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